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A  longstanding  problem  in  radar  and  electromagnetic  scattering  measure¬ 
ments  is  the  reconstruction  of  object  shape  and  detail,  i.e.  an  image,  from 
far  field  scatter  data  to  be  used  in  object  identification  and  classification. 
During  the  period  covered  by  this  report  we  have  been  able  to  demonstrate 
in  a  feasibility  study  the  first  reconstruction  of  a  3-D  image  of  a  perfectly 
reflecting  object  from  its  wave-vector  diversity  (multi frequency  and  multi¬ 
aspect)  scatter  data  making  use  of  a  new  Weighted  Fourier.  Domain  Projection  \ 
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'worm  (OJFVPT )  and  establish  a  relationship  between  holography  with  wavelength 
diversity  and  inverse  scattering.  Coherent  wave  vector  diversity  radar  tech¬ 
niques  are  used  to  access  a  finite  volume  of  the  3-D  Fourier  space  (also  known 
1  as  reciprocal  space  or  p-space)  of  the  scattering  object.  The  WFDPT  permits  the 
j.  use  of  hybrid  (digital/optical)  computing  that  enables  the  retrieval  and  display 
of  3-D  image  information  in  parallel  slices  or  cross-sectional  outlines.  The 
■  major  attributes  of  this  approach  as  compared  to  totally  digital  computing  are 
'  its  potential  for  displaying  a  true  3-D  image  in  real-time.  Wave-vector  di- 
I  versity  methods  appear  suitable  for  the  imaging  of  two  classes  of  practical  ob¬ 
jects  namely  non-dispersive  perfectly  reflecting  objects  of  the  type  often 
;  encountered  in  radar  (and  sonar)  and  semitransparent  weakly  scattering  objects 
such  as  certain  ultrasound  and  light  scattering  objects  encountered  in  biology 
“J  and  medic ine.^Jhe  method  has  several  unique  characteristics.  It  furnishes 
£  true  super- resolution,  i.e.  resolution  exceeding  the  classical  Rayliegh  Limit 
4  of  the  available  recording  aperture,  in  this  case  a  highly  thinned  (widely  dis- 
m  persed)  broad-band  coherent  receiver  array.  True  super-resolution  is  achieved 
because  of  an  inherent  aperture  synth'-sis  due  to  frequency  diversity  (frequency 
h  scanning,  stepping  or  comb  illumination)  and  conversion  of  spectral  degrees  of 
freedom  into  spatial  image  detail.  Accordingly,  when  applied  to  the  imaging  of 
ki  a  dispersive  object,  a  target  signature  rather  than  a  geometrical  image  should 
^  be  expected.  Such  a  target  signature  could  still  be  useful  in  object  identifi- 
cation  and  classification  since  it  contains  information  pertaining  to  the 
&  material  composition  of  the  object  intermixed  with  geometrical  image  detail. 
:-^?The  use  of  frequency  diversity  was  found  to  lead  to  a  unipolar  impulse  response. 

nThis  is  very  useful  in  suppressing  coherent  noise  (speckle)  which  is  known  to  be 
the  major  drawback  of  coherent  imaging. 

!  Preliminary  work  on  3-D  image  display  has  yielded  encouraging  results  on 

'3-D  display  from  a  series  of  weighted  projection  holograms  of  various  slices  of 
,  a  test  object.  The  projection  holograms  were  viewed  in  rapid  succession  using 
,  the  virtual  Fourier  transform. 

To  identify  optimal  and  practical  approaches  to  waye-vector  diversity 
data  acquisition  a  unique  microwave  measurement  system  has  been  assembled, 

’ v  .installed  and  tested  in  our  anechoic  chamber  facility.  The  system  was  used  also 
i  in  the  study  of  TDR  (target  derived  reference)  methods  in  which  a  reference  for 
phase  measurement  can  be  furnished  by  the  scattering  object  eliminating  thus  the 
r  heed  for  costly  local  oscillator  distribution  networks  and  eliminating  at  the 
-.same  time  undesirable  range  phase  ambiguities  from  the  collected  data.  The 
measurement  facility  is  computer  aided  furnishing  thereby  semi-automatic  control 
i .  of  object  positioning  or  orientation,  frequency  stepping,  data  acquisition  and 
[  storage,  and  final  data  correction  and  analysis.  A  high  resolution  CRTdisplay 
['*.  enables  the  display  of  weighted  projection  hologram  computed  from  the  p-space 
.  data  making  use  of  the  WFDPT.  Preliminary  results  of  this  measurement  system 
'capabilities  included  in  this  report  confirm  its  tremendous  versatility.  At 
this  stage,  the  program  strongly  suggests  the  practical  feasibility  of  a  new 
generation  of  cost  effective,  real-time,  super-resolving  3-D  imaging  radars  that 
can  because  of  their  3-D  image  slicing  characteristic  be  appropriately  referred 
\  to  as  Tomographic  Radon*  (Tomos=slice  in  Greek). 

"■  Finally,  a  study  of  3-D  imaging  using  other  forms  of  broadband  radiation 
”such  as  impulsive,  random  noise  and  particularly  thermal  emission  for  passive 
3-D  wavelength  diversity  imaging  has  been  also  initiated. 
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HIGH  RESOLUTION  FREQUENCY  SWEPT  IMAGING 


1.  Introduction 

The  aim  of  the  research  work  outlined  in  this  final  report  was  the  analysis 
and  investigation  of  methods  by  which  frequency  or  wavelength  diversity  techniques 
can  be  employed  to  impart  to  a  hicjhly  thinned,  and  therefore  cost-effective, 
longwave  (microwave  or  ultrasound)  imagina  aperture  resolution  capabilities  better 
than  its  monochromatic  classical  (Rayligh)  limit  achieving  thereby  super-resolution 
by  means  of  frequency  synthesized  apertures.  This  approach  to  longwave  imaging 
gains  practical  significance  when  one  considers  the  current  highly  developed  state 
of  the  art  of  broadband  microwave  gear  suitable  for  use  in  a  new  generation  of 
cost-effective  high  resolution  microwave  imaging  radars  utilizing  frequency 
diversity  techniques. 

It  is  well  known  that  the  development  of  longwave  holographic  imaging  systems 
possessing  resolution  and  image  quality  approaching  those  of  optical  systems  is 
hampered  by  three  factors:  (a)  prohibitive  cost  and  size  of  longwave  imaging 
apertures,  (b)  rapid  deterioration  of  longitudinal  resolution  with  range,  (c) 
inability  to  view  a  3-D  image  as  with  optical  Fresnel  holograms  because  of  a 
wavelength  scaling  problem  and  (d)  degradion  of  image  quality  by  speckle  or 
coherent  noise  because  of  the  low  numerical  apertures  attainable  with  present 
techniques.  For  example,  a  longwave  imaging  aperture  operating  at  a  wavelength 
of  3  cm  should  be  about  3  km  in  size  in  order  to  achieve  image  resolution  comparable 
to  an  ordinary  photographic  camera.  In  addition  to  inconvenient  size,  the  cost 
of  filling  such  a  large  aperture  with  suitable  coherent  sensors  is  clearly  prohibitive. 
Furthermore,  recall  that  in  conventional  longwave  holography  when  optical  image 
retrieval  is  utilized,  it  Is  necessary  to  store  the  longwave  hologram  data  (fringe 
pattern)  in  an  optical  transparency  suitable  for  processing  on  the  optical  bench 
using  laser  light.  In  order  to  avoid  longitudinal  distortion*  of  the  reconstructed 
image,  the  size  of  the  optical  hologram  replica  must  be  m  (=Xlong/Xlaser)  times 
smaller  than  the  longwave  recording  aperture.  For  the  example  cited  earlier,  this 
means  an  optical  hologram  replica  of  less  than  a  millimeter  in  size.  It  is 
certainly  not  possible  to  view  a  virtual  3-D  Image  through  such  a  minute  hologram 
even  with  optical  aids  since  these  tend  to  introduce  their  own  longitudinal 
distortion.  As  a  result,  longwave  holographers  have  long  learned  to  forgo  3-D 
Imagery  and  settled  Instead  for  2-D  Imagery  obtained  by  projecting  the  reconstructed 
real  Image  on  a  screen.  This  permits  lowering  of  the  reduction  factor  m  and 
consequently  relaxing  the  resolution  requirements  of  the  photographic  film  which 
allows  in  turn  the  use  of  highly  convenient  Polaroid  transparency  film  for 
preparation  of  the  optical  hologram  replica.  Because  of  the  small  size  (measured 
In  wavelength)  of  longwave  apertures  attainable  in  practice  and  the  above  methods 
of  viewing  the  real  image,  speckle  noise  is  always  present  leading  to  degration 
in  image  quality. 


*  Longitudinal  distortion  causes  for  example  the  image  of  a  sphere  to  appear 
elongated  in  the  range  direction  like  a  very  long  ellipsoid. 


In  this  report  we  summarize  the  main  results  of  our  investigation  under 
this  grant.  Our  findings  show  that  frequency  diversity  techniques  not  only 
circumvent  the  limitations  discussed  above  but  provide  a  means  of  viewing 
true  3-D  images  of  distant  objects  such  as  satellites  and  aircraft.  It  is  worth¬ 
while  to  point  out  that  our  studies  of  wave-vector  diversity  imaging  (or  frequency 
swept  imaging)  were  motivated  to  some  extent  by  evidence  of  super-resolved  "imaging" 
capabilities  in  the  dolphin  and  the  bat  which  are  known  to  use  frequency  swept 
(chirp)  signals  in  their  "sonar"  to  discern  small  objects  in  their  environment. 


2.  Summary  of  Important  Results 

The  main  findings  of  the  study,  details  of  which  are  given  in  the  appendices  and 
our  publications  (see  list  of  publications),  are  outlined  next. 

(a)  Wave-vector  diversity  (multi frequency  and  multiaspect)  techniques  can 
be  used  to  enhance  the  amount  of  object  information  collected  by  a  broadband 
coherent  aperture  deployed  in  the  far  field  of  the  scattering  object.  Thus  the 
data  collected  by  a  highly  thinned  array  of  coherent  receivers  intercepting  the 
wavefield  scattered  from  a  distant  3-D  reflecting  object,  as  the  frequency  of  its 
illumination  (and/or  its  direction  of  incidence)  are  changed  (see  Fig.  1-a  for 
example),  can  be  stored  as  a  3-D  data  manifold  in  p-space  (Fig.  2)  from  which  an 
image  of  the  object  can  be  retrieved  by  means  of  a  3-D  Fourier  Transform.  The 
size  and  shape  of  the  3-D  data  manifold,  and  therefore  the  resolution,  depend  on 
the  relative  positions  of  the  object, the  transmitter  (illuminator) ,and  the  receiving 
array  and  on  the  spectral  width  of  the  illumination  utilized. 

(b)  The  data  collected  must  be  corrected  for  a  quadratic  phase  factor  F 
(caused  by  the  unequal  distances  between  the  object  and  the  receiving  stations 
forming  the  widely  dispersed  imaging  array)  before  it  is  stored  in  a  3-D  manifold 
in  p-space  and  an  undistorted  image  of  the  3-D  reflecting  object  reconstructed 
through  the  3-D  Fourier  transform  operation.  A  bothersome  range-azimuth  ambiguity 
is  also  avoided  through  elimination  of  this  quadratic  phase  term. 

The  most  promising  methods  for  data  acquisition  and  correction  is  that  which 
utilizes  a  target  derived  reference  (TDR)  at  the  synchronous  detectors  of  the 
various  receivers  to  correct  for  the  unequal  phase  shifts  or  propagation  time 
delays  from  the  object  to  each  receiver.  In  this  approach  the  data  furnished  by 
the  various  receivers  of  the  recording  array  is  free  of  the  undesired  factor  F. 
Therefore  no  additional  processing  by  a  computer  will  be  necessary  before  filing 
the  data  at  the  appropriate  locations  in  p-space.  The  TDR  method  has  several 
advantages  which  include: 

(i)  Elimination  of  the  need  for  a  costly  and  unreliable  central 
local  oscillator  distribution  network. 

(ii)  Because  TDR  results  in  a  recording  configuration  similar  to 
that  of  a  lensless  Fourier  Transform  hologram,  the  resolution  requirements 
from  the  recording  device  are  greatly  relaxed*. 


♦A.  Macovski,  "Hologram  Information  Capacity",  J.  Opt.  Soc.  Am.,  Vol.  60, 
Jan.  1970,  pp.  21-29. 


In  longwave  holography  this  fact  is  translated  into  a  significant 
reduction  of  the  number  of  receiving  elements  in  the  recording  aperture. 
In  addition  the  use  of  TDR  allows  us  to  place  all  the  resolving  power 
of  the  recording  aperture  on  the  target.  This  means  that  high 
resolution  images  of  distant  isolated  targets  should  be  feasible  with 
array  apertures  consisting  of  tens  of  elements.  The  ability  to 
synthesize  a  2-D  receiving  aperture  with  a  Wells  arrayt  consisting  of 
two  orthogonal  linear  arrays  one  of  transmitters  and  the  other  of 
receivers  provides  further  means  of  reducing  the  number  of  stations 
needed  for  data  acquisition  without  sacrifice  in  resolution.  A 
frequency  swept  Wells  array  of  10  transmitters  and  10  receivers  using 
a  (2-4)  GHz  sweep  should  be  able  to  easily  furnish  104  3-D  distinguishabl 
resolution  cells  on  the  target  which  is  more  than  sufficient  for 
discerning  the  scattering  centers  on  practical  targets. 

( i v )  Greater  immunity  to  phase  fluctuations  arising  from 
turbulance  and  inhomogenieties  in  the  propagation  medium  because  both 
the  reference  and  imaging  signals  arriving  at  each  receiving  element 
of  the  aperture  travel  roughly  over  the  same  path. 

(v)  TDR  eliminates  the  range  azimuth  ambiguity  and  excessive 
bandwidth  problems  that  arise  in  fast  frequency  swept  imaging  when 

the  reference  signal  for  the  array  aperture  is  distributed  instead  from 
the  illumination  source  or  a  centrally  located  local  oscillator  phase 
locked  to  it. 

Two  TDR  methods  have  been  considered  to  some  extent  in  our 
work  to  date.  In  one  method  which  we  term  LFTDR  (Low  F K&que.ncy 
Target  VeAlve.d  Rej JeAence),  the  object  is  assumed  to  be  illuminated 
simultaneously  with  a  high  frequency  imaging  signal  and  a  low  frequency 
signal  that  is  a  subharmonic  of  the  illuminating  frequency.  The 
subharmonic  reference  frequency  u>  is  chosen  such  that  k  £«1,£  being 
the  maximum  linear  dimension  of  tne  object  and  k  =  w  /c,  c  being 
the  velocity  of  light.  This  places  scattering  from  tne  object  in 
the  Rayliegh  region  where  the  object  behaves  as  point  scatterer  with 
zero  phase  contribution.  The  far  field  phase  of  the  reference  signal 
at  any  receiver  is  therefore  entirely  due  to  propagation  between  a 
reference  point  formed  at  the  object  to  the  receiver.  A  method  for 
measuring  this  reference  signal  phase  and  using  it  to  correct  the 
imaging  signal  phase  due  to  propagation  has  been  proposed  by  Porter* 
and  analyzed  for  a  one-dimensional  object  geometry.  The  reference 
signal  phase  and  the  imaging  signal  phase  are  measured  separately 
at  each  receiving  station  with  the  aid  of  two  receivers  whose  local 
oscillators  (L.O)'s,  one  at  the  reference  frequency  and  one  at  the 
imaging  frequency,  are  phase-locked  only  to  each  other  and  not  to  a 
central  local  oscillator  as  would  be  the  case  were  we  to  use  a 
conventional  receiver  array.  Phase  locking  of  the  two  LO's  can  be 
accomplished  by  simply  making  the  imaging  L.O  a  harmonic  of  the 


+C.N.  Nilsen  and  D.N.  Swingler,  "Quasi-Real-Time  Inertialess  Microwave 
Holography",  Proc.  IEEE  (Letters),  Vol.  65,  March  1975,  pp.  491-492. 
*R.P.  Porter,  "A  Radar  Imaging  System  Using  the  Object  as  Reference", 
Proc.  IEEE  (Letters),  Vol.  59,  Feb.  1971,  pp.  307-308. 


reference  L.O.  This  would  eliminate  the  difficulties  encountered  in 
the  implementation  of  large  or  giant  thinned  coherent  receiving  arrays 
of  the  type  required  here,  namely  the  distribution  of  a  central  local 
oscillator  signal.  A  great  reduction  in  cost  and  effort  associated 
with  installation  of  a  central  L.O.  distribution  network  can  thus  be 
achieved.  This  cost  reduction  should  be  compared  however  with  the 
cost  of  implementing  a  LFTDR.  Because  of  the  large  difference  between 
the  high  frequency  imaging  frequencies  and  the  low  frequency  reference 
frequency  required  for  the  high  resolution  imaging  of  practical  objects, 
the  same  microwave  gear  can  not  be  used  for  both  frequencies.  This 
could  increase  system  cost.  In  addition  since  the  measured  reference 
phase  must  be  multiplied  by  a  factor  B  equal  to  the  ratio  of  the 
imaging  to  the  reference  frequency  before  being  used  as  a  reference 
phase  in  the  imaging  signal  measurement,  any  errors  in  the  reference 
phase  measurement  will  also  be  amplified  by  this  ratio.  The  precision 
of  the  reference  phase  measurement  and  phase  error  analysis  are 
important  and  will  have  therefore  to  be  examined  further. 


Another  TDR  methods  which  we  call  the  F Jizquwcy  DZiplac-zd.  TaAgvt 
VexLvzd  RetfeAence  (FDTDR)  also  shows  promise.  In  this  method,  the 
analytical  details  of  which  are  outlined  in  appendix  I,  the  object 
is  illuminated  simultaneously  during  the  sweep  with  two  phase  locked 
imaging  frequencies  w,  and  w~  ~  id,  +  Aw,  Aw  being  a  small  incremental 
frequency.  This  can  Be  realized  also  by  single  side  band  modulation 
of  the  swept  signal  or  by  phase  locking  two  sweep  oscillators. 
Measurement  of  the  differential  phase  between  the  signals  scattered 
from  the  target  at  these  frequencies  yields  Aw  /D  , 

“c  + 


f^)’  ^  being 


the  distance  from  the  transmitter  to  the  object  and  R  being  the 
distance  from  the  object  to  the  receiver.  Multiplication  of  this 
phase  by  w,/Aw  yields  the  phase  factor  F  at  frequency  w,  which  would 
be  used  to1correct  the  phase  measured  at  w-,.  At  first  look  this 
method  would  appear  to  still  require  a  reference  local  oscillator. 
This  however  is  not  so  since  the  procedure  outlined  above  need  not 
involve  explicit  phase  measurements  and  multiplications.  For 
example  by  mixing  the  two  received  signals  at  w,  and  w,  +  Aw  in  a 
square  law  detector  at  each  receiver  a  beat  signal  at  frequency  Aw 


is  derived  whose  phase  is  equal  to  Aw  + 

c  R 


* 


The  phase  shift  of 


this  signal  due  to  the  object  is  effectively  zero  because  the  wavelength 
at  Aw  is  much  larger  than  the  object  extent  making  it  behave  effectively 
as  a  point  scatterer.  Harmonic  mixing  of  the  signal  w.  received  at 
each  receiver  with  this  beat  signal  should  yield  the  corrected  p-space 
data  at  w. .  Because  of  the  small  difference  Aw  between  the  two 
frequencies  w.  and  in?  utilized,  the  effect  of  phase  errors  due  to  system 
and  atmospheric  propagation  could  be  more  completely  cancelled  in  this 
method  than  in  the  low  frequency  TDR  methods.  The  small  difference  Aw 
means  also  that  unlike  the  LFTDR  case  the  same  microwave  gear  (antennas, 
transmission  lines  and  other  microwave  circuit  components)  can  be  utilized 
in  the  handling  of  the  reference  and  imaging  signals.  A  variation  of 
the  TDR  technique  involving  double  side-band  modulation  is  also  possible 
and  appears  to  be  more  simple  to  implement  than  the  single 

side-band  method. 
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(c)  Because  in_addition  to  being  dependent  on  geometry , the  dimensions 
of  the  3-D  data  in  p-space  shown  in  Fig.  2  are  dependent  on  the  spectral 
range  of  the  illumination,  super-resolution  (i.e.  resolution  beyond  the 
classical  limit  of  the  available  physical  aperture)  is  achieved.  This 
aperture  synthesis  by  wave-vector  or  frequency  diversity  helps  cut  down 
array  cost  (since  a  thinned  array  can  be  used  to  frequency  synthesize  a 
large  array  with  higher  filling  factor). 

(d)  Foa'u.eA  domain  PKOjecXion  TheoaemA  (see  appendixes  II  and  III 
for  details)  enable  the  generation  of  two  dimensional  holograms  from 

projections  (or  weighted  projections)  of  the  corrected  p-space  3-D  data 
manifold  of  Fig.  2  permitting  thereby  optical  image  retrieval  of  the  3-D 
object  in  slices  parallel  to  the  projection  plane  one  at_a  time.  For 
example.  Fig.  1-b  shows  the  projection  hologram  for  the  p-space  data 
obtained  in  a  computer  simulation  of  the  arrangement  shown  in  Fig.  1-a. 

The  central  cross-sectional  outline  of  the  object  (the  two  1  m  diameter 
reflecting  spheres  of  Fig.  1-a)  retrieved  from  this  projection  hologram 
by  means  of  a  2-D  Fourier  transform  carried  out  on  the  optical  bench  is 
shown  in  Fig.  1-c.  A  similar  example  is  shown  in  Figs.  3  and  4.  Figure 
3  shows  a  second  test  object  consisting  of  3-D  distribution  of  a  set  of 
8  point  scatterers  with  locations  and  spacings  given  in  the  Figure. 

Figure  4  shows  the  projection  holograms  corresponding  to  the  three  slices 
of  the  object  containing  the  point  scatterers  and  the  image  retrieved 
from  each.  The  sweep  width  in  this  example,  as  in  the  previous  example, 
was  (2-4)  GHz  however  the  number  of  receivers  in  the  recording  array 
has  been  reduced  from  50  to  16.  These  computer  simulations  demonstrate 
that  a  3-D  (lateral  and  longitudinal)  resolution  of  the  order  of 
twenty  centimeters*is  easily  achieved  with  a  frequency  sweep  covering  only 
(2-4)  GHz  using  a  broad-band  array  of  16  receivers  and  one  transmitter. 
Wider-spectral  windows  should  yield  better  resolution.  It  is  worthwhile 
to  note  in  this  respect  that  commercial  microwave  sweepers  and  synthesizers 
are  available  with  a  spectral  coverage  of  (.1-25)  GHz  indicating  a 
potential  for  practical  resolutions  of  the  order  of  possibility  few 
centimeters  with  cosc-effective  broad-band  apertures  consisting  of  tens 
of  receivers  operating  with  one  central  illuminator. 

(e)  The  viewing  or  the  display  of  a  true  3-D  image  of  the  various 

slices  or  cross-sectional  outlines  should  be  possible  by  reconstruction 
of  the  various  projection  holograms  in  rapid  succession  while  projecting 
the  reconstructed  real  images  of  the  corresponding  slices  on  a  rapidly 
moving  projection  screen.  The  screen  would  be  displaced  rapidly  (together 
with  the  Fourier  transforming  lens)  on  the  optical  bench  in  the  axial 
directions  by  small  amounts  proportional  to  the  distances  between  the 
various  slices.  In  another  approach  we  have  found  that  the  2-D  vintual 
VoivueA  of  a  projection  hologram  can  be  carried  out  by  simply 

viewing  (with  the  unaided  eye)  a  transparency  containing  an  array  of 
reduced  replicas  of  the  projection  hologram  arranged  side-by-side  with  a 
point  source.  The  image  retrieved  in  this  fashion  would  lie  in  the 
plane  of  the  point  source.  This  approach  has  the  potential  for  3-D 
display  by  viewing  the  virtual  images  retrieved  from  a  series  of  projection 
holograms  corresponding  to  different  slices  or  cross-sectional  outlines 

*This  means  10^  distinguishable  3-D  resolution  cells  in  the  (2  x  2  x  2)m^ 
volume  of  the  assumed  object. 
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of  the  object  passed  in  front  of  the  eye  in  rapid  succession  while 
moving  the  reconstruction  point  source  axially  back  and  forth  at  a 
suitable  rate  of  incremental  axial  displacements.  A  proposed  electro- 
optical  scheme  that  permits  carrying  out  this  procedure  in  real-time 
using  a  rapidly  recyclable  spatial  light  modulator  (SLM)  operating  in 
a  reflection  mode  is  shown  in  Fig.  5.  The  computer,  the  high  resolution 
CRT  and  the  projection  optics  are  used  to  project  reduced  noncoherent 
images  of  the  various  projection  holograms  in  rapid  succession  on  the 
SLM  while  the  axial  position  of  the  reconstruction  point  sources  is 
altered  rapidly  also  under  computer  control.  The  point  source  need 
not  be  derived  from  a  laser  in  order  to  yield  an  image  but  could  also 
be  a  miniature  "grain  of  wheat"  light  bulb.  Details  of  this  task  are 
found  in  Appendix  III. 

(f)  As  seen  in  (e),  unlike  monochromatic  longwave  holographic 
imaging,  there  is  no  specific  scaling  requirements  imposed  on  the 
projection  holograms  in  order  to  avoid  longitudinal  distortion  in  the 
optical  reconstruction  circumventing  thus  the  wavelength  scaling 
problem. 


(g)  Because  of  the  broad  spectral  extent  of  the  illumination 
used  and  ability  to  display  the  reconstructed  image  in  separate  slices, 
speckle  or  coherent  noise,  which  is  known  to  plague  coherent  imaging 
systems,  is  suppressed  making  the  system  behave  in  as  far  as  image 
noise  is  concerned  like  a  noncoherent  imaging  system  but  at  the  same 
time  enjoy  the  superior  detection  characteristics  associated  with 
synchroneous  detection  techniques. 

(h)  The  broad-band  nature  of  the  imaging  process  also  helps 
suppress  undesirable  image  detail  that  could  arise  from  object  resonance o 
which  could  seriously  degrade  image  quality  in  a  monochromatic  imaging 
system. 


(i)  The  data  collected  at  every  receiver,  represents  after 
correction,  essentially  the  frequency  response  of  the  scattering  object 
measured  from  a  different  aspect  angle.  Assuming  the  scattering  process 
is  linear,  this  frequency  response  is  related  to  the  impulse  response  of 
the  object  by  a  Fourier  transform  (see  ref.  5  in  List  of  publications). 
This  suggests  that  impulse  illumination  can  be  utilized  instead  of 
freguency  swept  illumination.  When  this  is  done, the  3-D  data  manifold 
in  p-space  may  be  generated  by  Fourier  transforming  the  impulse  response 
at  each  receiver,  correcting  the  data  for  the  Factor  F  mentioned  in  (b), 
and  storing  the  result  in  the  appropriate  p-space  locations  for  each 
receiver.  The  resulting  p-space  volume  accessed  in  this  fashion  can  then 
be  employed  as  described  earlier  to  yield  3-D  image  information.  Impulse 
illumination  is  desirable  in  certain  instances  of  rapid  target  motion 
but  may  be  more  difficult  to  implement  than  frequency  swept  illumination. 
Since  the  impulse  reponse  of  a  time  invariant  linear  system  can  also  be 
deduced  from  white  noise  excitation  and  corellation  of  the  output 
response  with  the  input  as  described  elsewhere  in  more  detail  (see  5 
in  list  of  publications),  it  follows  that  the  techniques  described  in 
this  report  for  coherent  broadband  radiation  should  be  equally  applicable 
with  minor  signal  processing  modification  to  noise-like  broadband 
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radiation  including  passive  black-body  radiation. 


(j)  Experimental  verification  for  both  the  principle  of 
frequency  diversity  imaging  and  the  TDR  concept  were  obtained  with  the 
aid  of  a  semi -automated  network  analyzer  configured  and  installed  in  a 
recently  refurbished  anechoic  chamber  within  the  scope  of  this  program(  see 
Figs.  6  and  7).  This  versatile  system  is  capable  of  vector  (amplitude 
and  phase)  measurements  of  wavefields  scattered  from  test  objects 
situated  in  the  anechoic  chamber  over  any  frequency  range  lying  in 
the  ( . l-18)GHz  range  for  a  variety  of  polarizations.  A  test  object 
consisting  of  two  parallel  cylinders  25  cm  apart  each  5  cm  in  diameter 
and  50  cm  long  was  mounted  on  a  rotating  styrofoam  pedestal  that  is 
under  computer  control  and  illuminated  as  shown  in  Fig.  8.  The 
distance  from  the  center  of  the  object  to  the  illuminating  parabolic 
antenna  to  the  left  and  the  receiving  horn  feeding  the  network 
analyzer  was  2.5  m.  The  complex  frequency  response  of  this  object 
was  measured  in  the  (5-14)GHz  range  and  the  data  stored  for  128  object 
orientation  covering  360°.  The  stored  data  was  corrected  for  range- 
phase  with  a  synthetic  TDR  generated  in  the  computer  and  the  corrected 
data  displayed  and  photographed  yielding  the  frequency  swept  hologram 
shown  in  Fig.  9  (c).  The  image  retrieved  from  this  hologram  via  an 
optical  Fourier  transform  carried  out  on  the  optical  bench  is  shown  in 
Fig.  9  (d).  For  reasons  of  comparison  a  computer  simulation  of  this 
experiment  assuring  a  (2-18)GHz  sweep  was  performed.  The  resultant 
range-phase  corrected  hologram  and  the  image  retrieved  from  it  optically 
are  shown  in  Fin.  9  (a)  and  (b).  Further  detail  on  this  phase  of  the 
program  were  reported  in  an  MSc.  thesis  made  part  of  this  report  in 
Appendix  IV.  This  part  of  the  program  is  being  continued  with  the  aim 
of  further  enhancing  measurement  accuracy  and  demonstrating  imaging  of 
a  nonsimple  3-D  test  object  such  as  a  model  aircraft  utilizing  polarization 
diversity  to  further  enhance  image  quality. 


3.  Conclusions. 

The  primarily  analytical  and  numerical  study  of  frequency  diversity 
imaging  performed  under  this  grant  demonstrates  conclusively  the  feasibility 
of  a  new  generation  of  coherent  broadband  imaging  radars  capable  of 
furnishing  3-D  image  detail  of  distant  target  with  cost  effective  giant 
apertures  and  efficient  digital/optical  signal  processing. 

Future  work  in  this  area  will  focus  more  on  the  analysis  and  identification 
of  optimal  methods  for  data  acquisition,  processing  and  3-D  display.  The 
ultimate  aim  is  the  generation  of  design  criteria  for  a  prototype  system 
and  its  assessment  in  the  3-D  imaging  of  low  flying  aircraft  passing  within 
range  of  our  facilities  on  route  for  landing  at  the  Philadelphia  Airport. 
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Fig. 2  Three  dimensional  p-space  data  generated  by 
frequency  sweeping  and  collected  by  a  2-D 
array  of  receivers. 
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Projection  holograms  and  their  optical 
reconstructions  for  the  set  of  point  scatterers 
in  Fig. 7. 10  at  different  R’  planes.  (a)  Hologram 
and  reconstructed  image  of  scatterers  at  R*=-z 
plane,  (b)  Hologram  and  image  at  R,=0  plane. 

(c)  Hologram  and  image  at  R'=z  plane,  x  =y 
=zo=100cm.  20  °  ° 
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True  3-D  Image  reconstruction  based  on  the  virtual  Fourier  transform. 
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Block  Diagram  of  Automated  Microwave  Network  Analyzer  Showing  Interface  to 
MINC  II  Computer  via  IEEE  433  Standard  Interface.  Bun.. 


Microwave  Anechoic  chamber  shovjinc  illuminator  antenna 


Two  views  of  dual -cylinder  test  object  in  Anechoic  chamber, 
(a)  View  showing  illuminator  to  the  left  and  the  receiving 
horn  on  the  right  separated  by  absorbing  barrier,  (b)  View 
showing  test  object  mounted  on  rotating  styrofoam  pedastel. 
Cylinders  are5  cm  in  diameter,  50  cm  long,  25  cm  apart. 


Fig.  9.  Frequency  swept  holograms  and  retrieved  images  for  a  dual-cylinder 
test  object,  (a)  Computed  frequency  swept  hologram  for  a  (2-18)GHz 
sweep  and  (b)  retrieved  image;  (c)  measured  frequency  swept  hologram 
for  a  (5-14)GHz  sweep  and  (d)  retrieved  image. 
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APPENDIX  I 


The  Frequency  Displaced  Target  Derived  Reference 


A  second  TDR  method  which  we  refer  to  as  a  Vtie.qu.mcy  V-Upluced 
T&igeX  Vetu.ved  R evened  (FDTDR)  method  also  shows  promise.  This 
method  involves  simultaneous  illumination  of  the  object  with  two  phase 
locked  imaging  frequencies  to.  and  w?  =  to.  +  Ato  that  differ  by  a  small 
frequency  increment  Ato.  1  Referring  to  eq.  (9)  of  ref. 

10  (see  list  of  publications)  we  can  write  for  the  far  field  at  a  given 
receiver  location  RR, 


'|t^{k^  »RR) 


e'jk1<RT+RR)  /ufr) 


r 


dr 


(1) 


^2 ( ^2 ’ 


j(k1+Ak)  -jkj(RT+RR) 

2iRR—  6 


-jAk(RT+RR) 


<1+$r> 


(2) 


where  k^  2  =  Wj  2/c  and  Ak  =  Ato/c. 

By  making  Aw/ to  «1  the  integral  in  (2)  will  approach  that  in  (1).  The 
only  difference  between  the  far  fields  ip.  and  ip~  at  the  receivers  is  then 
the  phase  term  Ak(R^  +  RR).  Measurement  of  this  phase  difference  yields 
(R  +  Rr)  since  Ak  Is  known.  This  information  can  be  used  to  correct 
the  phase  of  either  the  tj>.  or  ip„  signals  to  obtain  the  required  p-space 
information. 


r(p)  =  /U(r)  e  j  p  *  r  dr 


(3) 
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Appendix  II 


HOLOGRAPHY,  WAVE-LENGTH  DIVERSITY 
AND  INVERSE  SCATTERING 


ABSTRACT 

The  use  of  wavelength  diversity  to  enhance  the  performance  of 
thinned  coherent  imaging  aperatures  is  discussed.  It  is  shown  that 
wavelength  diversity  lensless  Fourier  transform  recording  arrangements 
that  utilize  a  reference  point  source  in  the  vicinity  of  the  object 
can  be  used  to  access  the  three-dimensional  Fourier  space  of  non- 
dispersive  perfectly  reflecting  or  weakly  scattering  objects.  Hybrid 
(opto-digital )  computing  applied  to  the  acquired  3-D  Fourier  space 
data  is  shown  to  yield  tomographic  reconstruction  of  3-D  image  detail 
either  in  parallel  or  meridonal  (central)  slices.  •  Because  of  an 
inherent  ability  of  converting  spectral  degrees  of  freedom  into  spatial 
3-D  image  detail  true  super-resolution  is  achieved  together  with 
suppression  of  coherent  noise.  The  similarity  of  the  kev 
equations  derived  to  those  of  inverse  scattering  theory  is  pointed 
out  and  the  feasibility  of  using  other  forms  of  broadband  radiation 
such  as  impulsive,  noise  and  thermal  is  discussed.  Finally,  the 
potential  of  utilizing  wavelength  diversity  imaging  in  microscopy  and 
telescopy  are  discussed. 


INTRODUCTION 

A  frequently  encountered  question  in  the  science  of  image  formation 
is  how  to  make  an  available  aperature  collect  more  information  about 
the  scene  or  object  being  imaged  in  order  to  enhance  its  resolving 
power  beyond  the  classical  Rayleigh  limit.  This  process  is  known  as 
super- resolution  and  is  relevant  to  all  imaging  systems  whether  holo¬ 
graphic  or  conventional.  There  are  five  known  methods  for  achieving 
super-resolution.  These  include:  weighting  or  apodization  of  the 
1  2 

aperature  data  *  ;  analytic  continuation  of  the  wavefield  measured  over 

3  4  5 

the  aperature  *  ;  use  of  evanescent  wave  illumination  ;  maximum  entropy 

7 

method  ;  and  use  of  the  time  channel  .  Weighting  and  analytical 
continuation  techniques  are  known  to  become  rapidly  ineffective  as  the 
signal  to  noise  ratio  of  the  data  collected  decreases.  Maximum 
entropy  techniques  are  known  to  be  more  robust  as  far  as  noise  is 
concerned  but  involve  ususally  extensive  computation.  Illumination 
with  evanescent  waves  is  practical  in  limited  situations  where  full 
control  of  the  recording  arrangement  exists  as  in  microscopy  for  example. 


This  leaves  the  time  channel  approach  in  which  one  can  collect  in 
time  more  information  about  the  object  through  the  available  recording 
aperature  by  altering  the  object  aspect  relative  to  the  aperature  by 
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means  of  roation  or  linear  motion  *  or  by  altering  the  parameters  of 
the  illumation  such  as  directions  of  incidence,  wavelength  and/or 
polarization.  These  later  operations  are  known  to  increase  the  degrees 
of  freedom  of  the  wavefields  impinging  on  the  recording  aperature 
enhancing  thereby  their  ability  to  convey  information  about  the  nature 
of  the  scattering  object.  Sophisticated  imaging  systems  endevour  to 
convert  the  nonspatial  degrees  of  freedom  of  the  wavefield,  e.g., 
angular,  spectral  and  polarization  to  spatial  image  detail  enhancing 
thereby  the  resolution  capability  beyond  the  classical  Rayleigh  limit 
of  the  available  physical  aperature.  Obviously  such  procedures  involve 
more  signal  processing  than  that  performed  by  conventional  imaging 
with  lens  systems  or  holography. 

In  this  paper  we  consider  generalizing  the  holographic  concept  to 
include  wavelength  diversity  as  a  means  of  enhancing  resolution.  A 
quick  examination  of  the  basic  equations  of  holography  reveals  that 
the  lensless  Fourier  transform  hologram  recording  arrangement  is 
amenable  to  this  generalization.  This  conclusion  is  used  then  as  a 
starting  point  for  a  Fourier  optics  formulation  of  wavelength  diversity 
imaging  of  3-D  (three  dimensional)  nondispersive  objects.  The  results 
show  that  measurement  of  the  multiaspect  or  multistatic  frequency  (or 
wavelength)  response  of  the  3-D  object  permits  accessing  its  3-D 
Fourier  space.  The  resulting  formulas  are  identical  to  those  obtained 

from  a  multistatic  generalization  of  inverse  scattering10,11,12 
establishing  thus  a  clear  connection  between  holography  and  the  inverse 
scattering  imaging  problem.  The  inclusion  of  wavelength  diversity  in 
holography  is  shown  to  have  several  important  features:  (a)  the 
availability  of  the  3-D  Fourier  space  data  permits  3-D  image  retrieval 
tomographically  in  parallel  or  meridonal  (central)  slices  or  cross- 
sectional  outlines  by  the  application  of  Fourier  domain  projection 
theorems,  (b)  suppression  of  coherent  noise  and  speckle  in  the  retrieved 
image,  (c)  removal  of  several  longstanding  constraints  on  longwave 
(microwave  and  acoustical)  holography  such  as  the  impractically  high 
cost  of  the  aperatures  needed,  the  inability  to  view  a  true  3-D  image 
as  in  optical  holography  because  of  a  wavelength  scaling  problem,  and 
minimization  of  the  effects  of  resonances  on  the  object. 

WAVELENGTH  DIVERSITY 

We  start  by  inquiring  into  the  conditions  under  which  the  data 
from  N  holograms  of  the  same  nondispersi ve  object  recorded  over  the 
same  aperture,  each  at  a  different  wavelength,  can  be  combined  to 
yield  a  single  image  superior  in  quality  to  the  image  retrieved 
from  any  of  the  individual  holograms. 


One  approach  to  answering  the  question  posed  above  would  be  to 

determine  the  conditions  under  which  the  well  known  formulas13  for  the 
focusing  condition,  magnification  and  image  location  in  holography  can 
be  made  independent  of  wavelength.  This  quickly  leads  to  the  conclusion 
that  wavelength  independence  can  be  met  if  a  reference  point  source 
centered  on  the  object  is  used  and  proper  scaling  of  the  individual 
holograms  by  the  ratio  of  recording  to  the  reconstruction  wavelength 

is  performed  before  super-position15’24.  The  former  condition  is  that 

for  recording  a  lensless  Fourier  transform  hologram14  where  the 
presence  of  the  reference  point  source  in  the  object  plane  leads  to 
the  recording  of  a  Fraunhofer  diffraction  pattern  of  the  object  rather 
than  its  Fresnel  diffraction  pattern  because  of  the  elimination  of  a 
quadratic  phase  term  in  the  object  wavefield  in  the  recorded  hologram. 

This  is  known  to  result  in  a  highly  desirable  reduction  in  the 
resolution  required  from  the  hologram  recording  medium  and  is  therefore  of 
.  practical  importance  especially  in  nonoptical  holography.  More 
detail  of  the  processing  involved  in  combining  the  data  in  multi 

wavelength  hologram  can  be  found  elsewhere15. 

Additional  insight  into  the  process  of  attaining  super-resolution 
by  wavelength  diversity  is  obtained  by  considering  the  concept  of 

wavelength  or  frequency  synthesized  aperature  16~20.  The  synthesis 
of  a  one  dimensional  aperature  by  wavelength  diversity  is  based  on 
the  simple  fact  that  the  Fraunhofer  or  far  field  diffraction  pattern 
of  a  nondispersive  planar  object  changes  its  scale,  i.e.  it  "breathes", 
but  does  not  change  its  shape  (functional  dependence),  as  the  wave¬ 
length  is  changed.  A  stationary  array  of  broadband  sensors  capable 
of  measuring  the  complex  field  variations  deployed  in  this  breathing 
diffraction  pattern  at  suitably  chosen  locations  would  sense  different 
parts  of  the  diffraction  pattern  as  the  wavelength  is  altered 
collecting  thereby  mo^e  information  on  the  nature  of  the  diffraction 
pattern  and  therefore  on  the  object  that  gave  rise  to  it  than  if 
the  wavelength  was  fixed  (stationary  diffraction  pattern).  Each 
stationary  sensor  in  the  array  is  thus  able  to  collect  as  the  wave¬ 
length  is  changed,  and  the  breathing  diffraction  pattern  sweeps  over 
it,  the  same  set  of  data  or  information  collected  by  a  movable 
sensor  mechanically  scanned  over  the  appropriate  part  of  the  diffraction 
pattern  when  it  is  kept  stationary  by  fixing  the  wavelength.  Hence 
the  term  wavelength  or  frequency  synthesized  aperture. 

The  orientation  and  location  of  the  wavelength  synthesized  aperture 
for  any  planar  distribution  of  sensors  deployed  in  the  Fraunhofer 
diffraction  pattern  of  a  planar  object  and  the  retrieval  of  an  image 

from  the  data  collected  has  been  treated  earlier15’1^.  It  was  clear, 
however,  that  extension  of  the  wavelength  diversity  concept  to  the 
case  of  3-D  objects  is  necessary  before  its  generality  and  practical 
use  could  be  established. 
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For  this  purpose  we  considered  as  shown  in  Fig. 1(a)  an  isolated 
planar  object  of  finite  extent  with  reflectivity  D(pq),  where  pQ 

is  a  two  dimensional  position  vector  in  the  object  plane  (xQ»yo).  The 

object  is  illuminated  by  a  coherent  plane  wave  of  unit-amplitude  and 
of  wave  vector  R.  *=  k  I.  produced  for  example  by  a  distant  source 

located  at  Ry.  The  wavefield  scattered  by  the  object  is  monitored  at 
a  receiving  point  designated  by  position  vector  Rr  belonging  to  a 

recording  aperture  lying  in  the  far  field  region  of  the  object.  The 
receiving  point  will  henceforth  be  referred  to  as  the  receiver  and  the 

source  point  at  the  transmitter.  The  position  vectors  p0>  Ry  and  Rr 

are  measured  from  the  origin  of  a  cartesean  coordinate  system  (xQ,  yQ, 

zQ)  centered  in  the  object.  The  object  is  assumed  to  be  nondispersive 

i.e.,  D  is  independent  of  k.  However,  when  the  object  is  dispersive 
such  that  D(pQ,k)  =  D1(py)D2(k)  and  D2(k)  is  known,  the  analysis 

presented  here  can  easily  be  modified  to  account  for  such  object 
dispersion  by  correcting  the  data  collected  for  D2(k)  as  k  is  changed. 


Referring  to  Figure  1(a)  and  ignoring  polarization  effects,  the 
field  amplitude  at  Rr  caused  by  the  object  scattered  wavefield  may  be 

expressed  as. 


ik  “J*&i  •  ~ik  rR 

¥(k.RR)  =  g-  /D(p0)e  1  1  e _ _  dpQ 

rR 


(1) 


where  dpQ  is  an  abbreviation  for  dxQdyQ  and  the  integration  is  carried 
out  over  the  extent  of  the  object.  Noting  that  fy  «  pQ  -  Ry,  Ry  =  -Ryl^ 

2  _ 

and  using  the  usual  approximations  valid  here:  rR  ~  RR  +  PQ/2RR  -1R  .  pQ 
for  the  exponential  in  (1)  and  rR  RR  for  the  denominator  in  (1)  where 
IR  =  Rr/Rr  and  IR  =  Ic./k  are  unit  vectors  in  the  Rr  and  directions 
respectively,  one1can  write  eq.  (1)  as, 


ik  -J'k(RT  +  Rr)  -j  P  •  Pft  . 

*(k'RR>-MLe  R/D(P0)R  °«p0.  (2) 

K 

where  we  have  used  the  fact  that  the  observation  point  is  in  the  far 

p 

field  of  the  object  so  that  exp  ( - jk  p  /2RD)  under  the  integral  sign  can 

**  ^  ^  _  A 

be  replaced  by  unity.  In  eq.  (2),  p  =  k  (lk  -  IR)=  px  Ix  +  py  Iy  + 

Pz  Iz  is  a  three  dimensional  vector  whose  length  and  orientation  depend 


5 


Fig. 


Transmitter 


1. 
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central J  slice  and  cross  sectional  outline  c  shown. 
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on  the  wavenumber  k  and  the  angular  positions  of  the  transmitter  and 
the  receiver.  For  each  receiver  and/or  transmitter  present,  p  indicates 
the  position  vector  for  data  storage.  An  array  of  receivers  for 
example  would  yield  therefore  as  k  is  changed  (  frequency  diversity) 
or  as  K  (=klk  )  is  charged  (wave-vector  diversity)  a  3-D  data  manifold. 

The  projection  of  this  3-D  data  manifold  on  the  object  plane  yields 
ij,(k,RT)  because  P  •  P0  =  Pt  •  P0  =  Px*0  +  P^  where  px  =  k(Ik  -  IR)x 

and  p  =  k ( 1.  -  ID)  are  the  cartesian  components  of  the  projection 

Pt  of  p  on  the  object  plane.  Accordingly  eq.  (2)  can  be  expressed  as. 


*(k*  RR}  =  life 


e_jk(Rj  + 


Rr)/D(x 


yj 


‘j(pxxo 


Vc 


o’Jo 


dx  dy 
O  0 


(3) 


Because  of  the  finite  extent  of  the  object,  the  limits  on  the  integral 
can  be  extended  to  infinity  without  altering  the  result.  The  integral 
in  (3)  is  recognized  then  as  the  two  dimensional  Fourier  transform 

D(px,py)  of  D(xQ,yo).  It  is  seen  to' be  dependent  on  the  object  reflec¬ 
tivity  function,  the  angular  positions  of  the  transmitter  and  the 
receiver  and  on  the  values  assumed  by  the  wavenumber  k  but  is  entirely 
independent  of  range.  Information  about  D  can  thus  be  collected  by 
varying  these  parameters.  Note  that  the  range  information  is 
contained  solely  in  the  factor  F  =  jk  exp  [ - j k ( R j  +  RR)]/2nRR  preceeding 

the  integral.  The  field  observed  at  has  thus  been  separated  into 

two  terms  one  of  which,  the  integral  D,  contains  the  lateral  object 
information  and  the  other  F,  contains  the  range  information.  The 
presence  of  F  in  eq.  (3)  hinders  the  imaging  process  since  it  complicates 
data  acquisition  and  if  not  removed,  gives  rise  to  image  distortion 
because  RR  is  generally  not  the  same  for  all  receivers.  To  retrieve 
an  image  of  the  object  via  a  2-D  Fourier  transform  of  eq.  (3),  the 
factor  F  must  first  be  eliminated.  Holographic  recording  of  the 
complex  field  amplitude  given  in  (3)  using  a  reference  point  source 
located  at  the  center  of  the  object  will  result  in  the  elimination  of 
the  factor  F  and  the  recording  of  a  Fourier  transform  hologram.  This 

operation  yields  D  over  a  two  dimensional  region  in  the  p  ,p  plane. 

_  *  y 

The  size  of  this  region,  which  determines  the  resolution  of  the  retrieved 
image  depends  on  the  angular  positions  of  the  transmitter  and  the 
receiver  and  on  the  values  assumed  by  k,  i.e.  the  extent  of  the  spectral 
window  used.  The  later  dependence  on  k  implies  super-resolution  imaging 
capability  because  of  the  frequency  synthesized  dimension  of  the  2-D 
data  manifold  generated.  Because  of  the  dependence  of  resolution  on 
the  relative  positions  of  the  object,  the  transmitter,  and  receiving 
aperture,  the  impulse  response  is  clearly  spatially  variant.  In  fact 
a  receiver  point  situated  at  for  which  p  is  normal  to  the  object 
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plane  can  not  collect  any  lateral  object  information  because  for  this 
condition  (p  .  pQ  =  0)  the  integrals  in  (2)  and  (3)  yield  a  constant. 

Such  receiving  point  is  located  in  the  direction  of  specular  reflection 
from  the  object  where  the  diffraction  pattern  is  stationary  i.e.  does 
not  change  with  k.  In  this  case  the  observed  field  is  solely  propor¬ 
tional  to  F  containing  thus  range  information  only.  Obviously  this 
case  can  easily  be  avoided  through  the  use  of  more  than  one  receiver 
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which  is  required  anyway  when  2-D  or  3-D  object  resolution  is  sought  *  . 

The  analysis  presented  above  can  be  extended  to  three  dimensional 

objects  by  viewing  a  3-D  object  as  a  collection  of  thin  merdional  or 

central  slices  as  depicted  in  Fig.  1(b)  each  of  which  representing 

a  two  dimensional  object  of  the  type  analyzed  above.  With  the  n-th 

slice  we  associate  a  cartesean  coordinate  system  x  ,y  ,z  that 

o  ■'o^  o^ 
n  n 

differ  from  other  slices  by  rotation  about  the  common  x  axis.  Since 

o 

the  vectors  p,  R^  and  RR  are  the  same  in  all  n-coordinate  systems, 

eq.  (3)  holds,  ip  (k,R  )  is  then  obtained  from  projection  of  the  three 
n  k 

dimensional  data  manifold  collected  for  the  3-D  object  on  the  xo,yQn 

plane  associated  with  the  n-th  slice.  An  image  for  each  slice  can 
then  be  obtained  as  described  before.  An  inherent  assumption  in  this 
argument  is  that  all  slices  are  illuminated  by  the  same  plane  wave. 

This  is  a  reasonable  approximation  when  the  3-D  object  is  weakly 
scattering  and  the  Born  approximation  is  applicable  or  when  the  3-D 
object  is  perfectly  reflecting  and  does  not  give  rise  to  multiple 
reflections  between  its  parts.  In  the  later  case  the  two  dimensional 
meridonal  slices  Dn(pQ  )  deteriorate  into  contours,  such  as  C  in 

Fig.  1(b)  defined  by  tRe  intersection  of  the  meridonal  planes  with 
the  illuminated  portion  of  the  surface  of  the  object.  Accordingly 
we  can  write  for  the  n-th  meridonal  slice  or  contour. 


,  -  ~JP • P_ 

*n(k.RR)  =  F/Dn(po  }  6  n  dpo  (4) 

n  n 

We  can  regard  D  (5  )  as  the  n-th  meridonal  slice  or  contour  of  a 

n  o 

three  dimensional  oBject  of  reflectivity  U(r)  where  r  is  a  three 
dimensional  position  vector  in  object  space.  This  means  that 
Dn(p0  )  =  U(r)  6(zq  )  where  6  is  the  Dirac  delta  "function". 

n  n 

Consequently  eq.  (4)  becomes, 
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~JP  •  p 

«i»n(k,RR)  =  F  /  U(r)  6  (*  ■)  e  °n  dp 

...  ,  .  “JP-?  -  (5) 

=  F  /  U(r)  6  (zQ  )  e  dr 

n 

where  dr  designated  an  element  of  volume  in  object  space  and  where 
the  last  equation  is  obtained  by  virture  of  the  sifting  property 
of  the  delta  function. 

Summing  up  the  data  from  all  slices  or  contours  of  the  object 
we  obtain. 


-jp.r 

Z  4)  =  F  /  U(r)  e  dr  =  tj/(p) 

n  n 


(6) 


because 


\ 


£  U(r)  6  (z  )  =  U(r). 
n 


Assuming  that  the  Factor  F  in  eq.  (6)  is  eliminated  as  before, 
equation  (6)  reduces  to 


-jp.r 

i|4p)  =  /  U(r)  e  dr 


(7) 


which  is  the  3-D  Fourier  transform  of  the  object  reflectivity  U(r). 
Wavelength  diversity  permits  therefore  accessing  the  3-D  Fourier  space 
of  a  nondispersive  object  providing  thereby  the  basis  for  3-D  Lensless 
Fourier  transform  holography.  An  alternate  formulation  to  that  given 
above  of  super-resolved  wave-vector  diversity  imaging  of  3-D  perfectly 
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conducting  objects  is  possible  by  extending  the  formulation  of  the 

inverse  scattering  imaging  problem*®’**  to  the  multistatic  case,  along 

12 

lines  that  are  similar  but  somewhat  different  than  those  given  by  Raz  . 
The  resulting  scalarized  formulas  are  identical  to  (7)  establishing 
thus  the  connection  between  the  holographic  and  the  inverse  scattering 
approaches  to  the  imaging  problem. 


THREE  DIMENSIONAL  IMAGE  RETRIEVAL 

The  above  considerations  of  multiwavelength  holography  have  lead 
us  to  determining  a  means  by  which  the  3-D  Fourier  space  of  the  object 
can  be  accessed  employing  synchroneous  detection.  It  is  clear  that 
once  the  3-D  Fourier  space  data  is  available,  3-D  image  detail  can  be 
retrieved  by  means  of  an  inverse  3-D  Fourier  transform  .which  can  be 
carried  out  digitally.  Alternately,  holographic  techniques 


<-xo.fc,*o> 


3.  Projection  holograms  and  their  optical  reconstructions  for 
the  set  of  point  scatterers  in  Fig.  2  at  different  planes 

(a)  Hologram  and  reconstructed  image  of  scatterers  at  R^-Zq 

plane,  (b)  Hologram  and  image  at  R^=0  plane,  (c)  Hologram 

and  image  at  R'~z„  plane,  x  =y„-z  =100  cm. 
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Fig.  4.  Arrangement  used  In  computer  simulation  of  wavelength  diversity 
tmaaina .  * 
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can  be  invoked  again.  Fourier  domain  projection  theorems  that  are 

dual  to  the  spatial  domain  projection  theorem^’2**  can  be  applied  to 
the  Fourier  space  data  to  produce  a  series  of  projection  holograms 
from  which  2-D  images  of  meridonal  or  parallel  slices  of  the  object 

can  be  retrieved  on  the  optical  bench  .  This  procedure  does  not 
involve  any  specific  scaling  of  the  size  of  the  optical  hologram 
transparency  relative  to  the  size  of  the  original  recording  aperture 
by  the  ratio  of  the  recording  to  the  reconstruction  wavelengths  as 
in  longwave  holography  where  the  scaling  necessary  for  viewing  a  3-D 
image  free  of  longitudinal  distortion  ususally  leads  to  an  impractically 
minute  equivalent  hologram  transparancy  that  cannot  be  readily  viewed 
by  an  observer.  The  lateral  and  longitudinal  resolutions  in  the 
retrieved  image  depend  now  on  the  dimensions  of  the  volume  in  Fourier 
space  accessed  by  wavelength  diversity.  This  volume  depends  on  the 
wavelength  range  and  on  the  recording  geometry.  Thus  the  longitudinal 
resolution  does  not  deteriorate  now  as  rapidly  with  range  as  in 
conventional  monochromatic  imaging  systems. 

/ 

An  example  of  computer  simulations  of  frequency  diversity  holo¬ 
graphic  imaging  of  a  3-D  object  consisting  of  eight  point  scatterers 
distributed  as  shown  in  Fig.  2  is  given  in  Fig.  3.  Shown  in  Fig.  3  are 
three  weighted  Fourier  domain  projection  holograms  and  the  corresponding 
optically  retrieved  images  for  three  equally  spaced  parallel  slices 
of  the  object  containing  distinguishable  2-D  distributions  of  scatterers. 
The  simulated  recording  arrangement  shown  in  Fig.  4  consisted  of  an 
array  of  16  receivers  equally  distributed  on  an  arc  extending  from  <j>  = 

40°  to  <{>  =  77.5°  surrounding  a  central  transmitter  capable  of  providing 
plane  wave  illumination  of  the  object.  The  results  shown  were  obtained 
with  microwave  imaging  in  mind  assuming  a  frequency  sweep  of  (2-4)GHz. 
They  clearly  indicate  a  leteral  and  longitudinal  resolution  capability 
of  the  order  of  25  cm.  Wider  sweep  widths  yield  better  resolution. 

For  example  a  (l-18)GHz  sweep  would  yield  a  3-D  resolution  of  the  order 
of  1.5  cm. 


DISCUSSION  AND  CONCLUSIONS 

Seeking  means  by  which  the  information  content  in  a  hologram  can 
be  increased  for  example  by  wavelength  diversity  we  have  arrived  at  a 
formulation  of  3-D  Lensless  Fourier  transform  holography  capable  of 
furnishing  3-D  image  detail  tomographically.  This  ability  of 
producing  3-D  images  in  slices  from  coherently  detected  wavefields 
enable  us  to  regard  the  method  also  as  coherent  tomography.  The  Fourier 
space  accessed  in  the  above  fashion  by  wavelength  diversity  can  be 
viewed  as  a  generalized  3-D  hologram  in  which  one  dimension  has  been 
synthesized  by  wavelength  diversity.  Such  a  generalized  hologram 
contains  not  only  spatial  amplitude  and  phase  data  as  in  conventional 
holography  but  also  spectral  information  and  hence  can  yield  better 
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resolution  than  the  classical  Rayleigh  limit  of  the  available  aperture 
operating  at  the  shortest  wavelength  of  the  spectral  window  used.  This 
super-resolving  property  is  further  enhanced  through  an  inherent  sup¬ 
pression  of  the  effects  of  object  resonances  and  coherent  noise  in  the 
retrieved  image,  the  latter  being  so  because  frequency  diversity  tends 
to  make  the  impulse  response  of  the  system  unipolar  resembling  that  of 
a  non-coherent  imaging  system  that  is  free  of  speckle  and  coherent  noise 
15 

artifacts  .  Futher  enhancement  of  information  content  and  resolution 
can  be  achieved  by  polarization  diversity  where  the  p  space  can  be 
multiply  accessed  for  different  nonredundant  polarizations  of  the 
illumination  and  the  receivers  and  the  resulting  polarization  diversity 
images  added  either  coherently  or  non-coherently  in  order  to  achieve  a 

15 

degree  of  noise  averaging  as  discussed  elsewhere  . 

The  removal  of  several  longstanding  constraints  on  conventional 
longwave  (microwave  and  acoustic)  holography  attained  through  the  use 
of  wavelength  diversity  as  described  here  leads  to  a  new  class  of  imaging  systems 
capable  of  converting  spectral  degrees  of  freedom  into  3-D  spatial  image 
detail  furnishing  thereby  true  super-resolution.  Wavelength  diversity 
is  applicable  to  the  imaging  of  two  classes  of  objects:  perfectly 
reflecting  objects  of  the  type  encountered  in  radar  and  sonar  and  weakly 
scattering  objects  of  low  or  known  dispersion  of  the  type  encountered 
in  biology  and  medicine.  The  practical  application  of  the  concepts 
presented  here  to  optical  wavefield  is  presently  under  consideration. 

The  availability  of  tunable  dye  lasers  and  electronic  imaging  devices 
suggest  interesting  possibilities  of  three  dimensional  wavelength 
diversity  microscopy.  Here  one  can  conceive  of  an  arrangement  in 
which  a  minute  semitransparent  object  with  homogeneous  or  known  dispersion 
is  transilluminated  by  a  collimated  coherent  light  beam  from  a  tunable 
dye  laser  which  can  also  be  made  to  provide  a  coherent  reference  point 
source  in  the  immediate  vicinity  of  the  object.  The  resulting  reference 
and  the  object  scattered  wavefields  are  intercepted  by  the  photocathode 
of  an  electronic  imaging  device  of  known  spectral  response  such  as  a 
vidicon.  Because  of  the  minute  size  of  the  object,  the  photocathode 
can  easily  be  situated  in  the  far  field  of  the  object.  Thus  nearly  a 
lensless  Fourier  transform  hologram  recording  arrangement  results.  The 
spatial  frequency  content  in  the  resulting  hologram  is  therefore  expected 
to  be  sufficiently  low  to  be  resolved  by  a  high  resolution  electronic 
imaging  device.  By  recording  and  digitally  storing  the  resulting  detected 
hologram  fringe  pattern  as  a  function  of  dye  laser  wavelength  onecan 
gain  access  to  the  3-D  Fourier  space  of  the  object  since  1^  and 

for  the  recording  geometry  are  precisely  known. 

A  similar  recording  arrangement  can  be  envisioned  in  the  active 
coherent  imaging  of  a  distant  reflecting  object  (active  telescopy)  where 
the  object  can  be  made  to  furnish  a  reference  point  source  situated  on 
its  surface  like. a  wavelength  independent  stationary  glint  point  or 
an  intentionally  placed  retroreflector.  Because  in  such  an  arrangement 
the  reference  and  the  object  wavefields  travel  over  the  same  path, 
atmospheric  effects  are  expected  to  be  minimized.  The  generation  of  an 
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object  derived  reference  geometry  in  longwave  (microwave  and  acoustic) 

20  27 

wavelength  diversity  imaging  has  been  described  elsewhere  5 

Finally  it  is  worthwhile  to  note  that  since  the  scattering  process 
is  linear  the  multiaspect  or  multistatic  frequency  or  wavelength 
response  measurements  referred  to  in  this  paper  can  be  obtained  also 
by  measuring  the  multi  aspect  impulse  response  followed  by  Fourier 

19 

transformation  of  the  individual  impulse  responses  measured  .  This 
means  that  impulsive  illumination  can  also  be  utilized.  Because  the 
impulse  response  of  a  linear  system  can  be  measured  by  using  random 

19 

noise  excitation  and  cross-correlating  the  output  with  the  input  ,  a 
possibility  of  using  random  noise  (white  light)  illumination  and 
cross-correlation  detection  techniques  as  a  means  for  accessing  the 
3-D  Fourier  space  of  the  object  also  emerges. 
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APPENDIX  III 


THE  VIRTUAL  FOURIER  TRANSFORM 
AND  ITS  APPLICATION  IN  THREE  DIMENSIONAL  DISPLAY 


ABSTRACT 


In  contrast  to  the  well  known  and  widely  used  instantaneous 
Fourier  transforming  property  of  the  convergent  lens  in  coherent 
(laser)  light,  the  "Virtual  Fourier  Transform"  (VFT)  capability  of 
the  divergent  lens  is  less  widely  known  or  used  despite  many  advantages. 
We  will  review  the  principle  of  the  VFT  and  discuss  its  advantages  in 
certain  applications.  In  particular  a  method  for  viewing  the  virtual 
Fourier  transform  of  a  two  dimensional  function  with  the  naked  eye 
using  an  ordinary  point  source  will  be  presented.  A  scheme  for  three- 
dimensional  image  display  based  on  a  "Fourier  domain  projection  theorem" 
utilizing  varifocal  VFT  is  described  and  a  discussion  of  the  properties 
of  the  displayed  image  given. 


INTRODUCTION 


Several  sophisticated  three  dimensional  (3-D)  imaging  techniques 

1  2  2 
such  as  x-ray  tomography  ,  electron  microscopy  ,  crystallography  , 

3 

wave-vector  diversity  imaging  and  inverse  scattering  ,  involve  measure¬ 
ments  that  give  access  to  a  finite  volume  in  the  3-D  Fourier  space  of 
a  3-D  object  function.  A  3-D  image  of  the  original  object  can  then 
be  reconstructed  by  computing  the  inverse  3-D  Fourier  transform.  The 
retrieved  iamge  normally  represents  the  spatial  distribution  of  a 
relevant  parameter  of  the  object  such  as  absorption,  reflectivity, 
scattering  potential,  etc. 

Obviously,  the  required  inverse  transform  can  be  performed  digitally. 
Digital  techniques  however  often  preclude  real-time  operation  partic¬ 
ularly  when  the  object  being  imaged  is  not  simple  but  contains  consider¬ 
able  resolvable  intricate  detail.  More  importantly,  because  of  the 
inherent  two  dimensionality  of  CRT  computer  displays, direct  true  3-D 
image  display  is  not  possible.  Present  day  computer  graphic  displays 
are  capable  of  displaying  3-D  image  detail  either  in  seperate  cross- 
sections  or  slices,  or  in  a  computed  perspective  (isometric)  view  of 
the  object,  or  in  some  instances  stereoscopically  where  an  illusion  of 
a  3-D  scene  is  created  in  the  mind  of  the  observer  who  is  required 
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usually  to  use  special  viewing  glasses  ’  . 


2 


Hybrid  (opto-digital )  computing  techniques  offer  an  alternate 
approach  to  3-D  image  retrieval  from  3-D  Fourier  space  data.  They 
furnish  as  shown  in  this  paper  the  ability  to  display  true  3-D  image 

2  3 

detail.  The  approach  is  based  on  "Fourier  Domain  Projection  Theorems"  * 
that  are  dual  to  "Spatial  or  Object  Domain  Projection  Theorems"  used  in 

radio-astronomy’^and  tomography*.  These  theorems  permit  the  recon¬ 
struction  of  3-D  image  detail  tomographically*  i.e.  in  slices  from 

2  3 

2-D  projections  of  the  3-D  Fourier  space  data  ’  .  Although  the  required 
2-D  Fourier  transform  can  be  carried  out  digitally,  the  emphasis  in 
this  paper  is  on  coherent  optical  techniques  for  performing  the  2-D 
Fourier  transform. with  particular  attention  to  implementations  that 
permit  the  execution  of  the  necessary  2-D  optical  transforms  of  the 
various  projection  hologram  sequentially  in  real-time.  Specific 
attention  is  given  to  a  technique  that  Utilizes  the  virtual  Fourier 
transform  which  permits  the  viewing  of  a  virtual  3-D  image  in  real-time. 

FOURIER  DOMAIN  PROJECTION  THEOREMS 

There  are  two  Fourier  domain  projection  theorems.  One  leads  to 
tomographic  object  reconstruction  in  parallel  slcies  and  is  called  the 
"weighted  Fourier  domain  projection  theorem"  the  other  leads  to 
tomographic  object  reconstruction  in  meridonal  or  central  slices  and 
can  therefore  be  called  the  "meridonal  or  central  slice  Fourier  domain 
projection  theorem" 

We  begin  by  considering  a  3-D  object  function  f(r)  with  r  =_ 
xlx  +  yly  +  zl2  being  a  position  vector  in  object  space.  Let  F(w) 

be  the  3-D  Fourier  transform  of  f(r)  defined  by, 

-jw.r 

F(w)  =  /  f(r)  e  dr  (1) 


where  dr  =  dx  dy  dz  and  w  =  w  I  +  w  I  +  w  I  is  a  position  vector 

x  x  y  y  z  z 

in  the  Fourier  or  spatial  frequency  domain. 

Consider  next  the  projection  of  F(w)  on  the  wY,  w„  plane  defined  by, 

a  jr 

Fp(WX’Wy)  =  (,  F{^  dWZ  • 

and  combining  eq.  (1)  and  (2), 

-j(w  x  +w.y  +w  z) 

Vwx’V  V  f<x>y>z)e  r  dxdydz}dwz  (3) 

★  '  , 

From  the  Greek  work  Tomos  meaning  slice. 
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Integrating  with  respect  to  wz  first  and  assuming  that  the  volume  in 
w  space  occupied  by  F(w)  is  sufficiently  large  we  obtain, 

-j(w  x  +  w  v) 

VWx’V  =  f^x,y,z^  8  e  ^  dxd^dz  (4) 

- j (w  X  +  w  v) 

-  {  f(x.y.o)  e  r  dxdy  (5) 

The  2-D  Fourier  domain  projection  Fp(wx»wy)  and  the  central  slice 

f(x,y,o)  through  the  object  form  thus  a  Fourier  transform  pair.  This 
may  be  symbolically  expressed  as, 

Fp(wx,Wy)  ++  f(x,y,o)  (6) 

Other  parallel  slices  through  the  object  at  z  =  zn,  zR  being  a 

constant  describing  the  z  coordinate  of  the  n-th  parallel  slice,  can 
in  a  similar  manner  be  related  to  "weighted"  Fourier  domain  projections 
of  F(w)  defined  by. 


Making  use  of  eq.  (1)  and  again  performing  the  integration  with  respect 
to  wz  first  we  obtain. 


(wx>wy)'e>f(x.y»zn) 


(8) 


which  indicates  that  the  weighted  projection  F  (w  ,w  )  and  the  n-th 

x  y 

parallel  object  slice  f(x,y,zn)  form  a  Fourier  transform  pair. 

Equation  (6)  is  seen  to  be  a  special  case  of  eq.  (8)  when  zn  *  o. 

Given  the  3-D  Fourier  space  data  manifold  F(w)  one  can  digitally 

compute  and  display  a  set  of  "weighted  projection  holograms"  F  (w  ,w  ). 

P » n  x  y 

A  corresponding  set  of  images  of  parallel  slices  or  cross-sectional 
outlines  of  the  3-D  object  can  then  be  retrieved  via  2-D  Fourier 
transform  operations  which  can  most  conveniently  be  carried  out 
optically  from  photographic  transparency  records  of  the  weighted 
projection  holograms  displayed  by  the  computer. 


i 
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Returning  to  eqs.  (1)  and  (2)  one  can  also  show  that  projections 

of  F(w)  on  arbitrarily  oriented  planes  other  than  the  w  ,w  plane 

x  y 

chosen  for  eq.  (2),  yields  "meridonal  projection  holograms"  that  are 
2-D  Fourier  transforms  of  corresponding  merdional  (central)  slices 
of  the  object.  This  is  the  "meridonal  Fourier  domain  projection 
theorem.  It  furnishes  the  basis  for  angular  multiplexing  of  the 
resulting  meridonal  projection  holograms  into  a  single  composite 
hologram  which  can  be  used  to  form  a  3-D  iamge  of  the  object  in  a 

Q 

manner  similar  to  that  in  integral  holography  which  is  increasingly 
being  referred  to  as  Cross  holography*. 

THE  VIRTUAL  FOURIER  TRANSFORM 


In  contrast  to  the  well  known  spatial  Fourier  transforming  pro¬ 
perty  of  the  convergent  lens  widely  used  in  coherent  optical  computing, 
the  complementary  virtual  Fourier  transform  capability  of  a  divergent 

lens10  is  less  widely  known  or  used  despite  many  attractive  features. 
This  is  surprising  since  the  power  spectrum  associated  with  the  VFT 
is  a  phenomenon  that  is  frequently  observed  in  daily  life  when  one 
happens  to  look  at  a  distant  point  source  such  as  a  street  light 
through  a  fine  mesh  screen  or  the  fine  fabric  of  transparent  curtain 
material.  The  spectrum  of  the  screen  transmittance  appears  then  as 
a  virtual  image  in  the  plane  of  the  point  source. 


The  VFT  concept  of  the  divergent  lens  Is  easily  derived  from  the 
Fourier  transform  expression  of  the  convergent  lens.  Figure  1  illustrates 
the  well  known  process  of  forming  a  real  Fourier  transform  with  a 
convergent  lens.  The  object  transparency,  with  complex  transmittance 
t(x  ,y  ),  is  placed  at  a  distance  d  in  front  of  a  convergent  lens  of 
focal  length  F  and  illuminated  with  a  normally  incidnet  collimated 
laser  beam.  The  complex  field  amplitude  of  the  wavefield  in  the 
back  focal  plane,  the  transform  plane,  is  given  by  the  well  known  formula 


T(x,y) 


J.  -jlF  1(1  "  F)(x2  +  y2)l 
If  e 


x  //  t(xQ,y0)  e 


•fr  <*  *„  ♦  y  V 


dx  dy^ 
o  o 


(1) 


in  which  the  integral  is  recognized  as  the  two  dimensional  Fourier 
transform  of  the  object  transmittance.  T(x,y)  becomes  the  exact 
Fourier  transform  of  t(xQ,yo)  when  d  =  F  that  is  when  the  object 

transparency  is  placed  in  the  Front  focal  plane  of  the  lens.  The 
power  spectrum  associated  with  the  transform  is  real  and  can  be 
projected  on  a  screen  placed  in  the  back  focal  plane.  It  is  also 
well  known  that  a  scaled  version  of  the  transform  can  be  obtained 
in  the  back  focal  plane  by  placing  the  object  tranparency  in  the 

g 

converging  laser  beam  to  the  right  of  the  lens  . 

*Named  after  Lloyd  Cross  the  originator  of  integral  holography. 


A 


r 


Object  plane 

y„  k 


Real  Fourier  transform 
w  Plane 


Collimated 
Laser  beam 


x 

/ 


h~ 


Fig.  1.  Real  Fourier  transform  formed  with  a  convergent  lens 
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Fig.  2.  Virtual  Fourier  transform  formed  with  a  divergent  lens 


Noting  that  eq.  (1)  does  not  change  when  we  replace  d  by  -d, 

F  by  -F,  xQ  and  yQ  by  -xQ  and  -yQ  respectively,  we.  can  arrive  at  the 

complementary  VFT  arrangement  illustrated  in  Fig.  2.  An  inverted 
transparency  t(xQ,yo)  is  placed  now  in  the  divergent  coherent  beam 

to  the  right  of  the  divergent  lens  {of  focal  length  -F)  and  a  VFT 
given  by  eq.  (1)  is  observed  in  the  virtual  focal  plane  of  the  lens. 

The  same  VFT  can  be  seen  by  removing  the  divergent  lens  and  replacing 
the  laser  beam  with  a  point  source  placed  at  the  origin  of  the  VFT 
plane  as  depicted  in  Fig.  3.  Thus  a  simple  way  of  viewing  the  power 
spectrum  associated  with  the  VFT  of  a  given  diffracting  screen 
(which  is  usually  a  Fourier  transform  hologram  or  a  projection 
hologram  of  the  type  described  above)  is  to  hold  the  screen  close 
to  the  eye  and  look  through  it  at  a  distant  bright  point  source.  The 
point  source  used  need  not  be  derived  from  a  laser.  In  fact  it  is 
preferable  for  safety  purposes  to  use  an  LED  or  a  spectrally  filtered 
minute  white  light  source  such  as  a  "grain-of-wheat"  subminiature 
incandescent  lamp  or  a  miniature  Christmas  tree  decorating  lamp 
covered  by  a  color  or  interference  filter.  This  has  the  added 
advantage  of  furnishing  a  measure  of  control  over  the  coherence  proper¬ 
ties  of  the  wavefield  impinging  on  the  screen  providing  thereby  a 
means  for  reducing  coherent  noise  in  the  observed  VFT  and  also,  as 
will  be  discussed  below,  a  means  for  coherent  or  noncoherent  super¬ 
position  of  VFT's.  As  the  distance  of  the  point  source  from  the 
diffracting  screen  is  decreased  in  order  to  make  it  compatible  with 
typical  laboratory  or  optical  bench  dimensions,  the  size  of  the 
observed  VFT  decreased  because  of  the  change  in  the  curvature  of  the 
wavefield  illuminating  the  diffracting  screen.  To  compensate  for  this 
effect  it  is  necessary  to  reduce  the  size  of  the  diffracting  screen  or 
transparency  often  to  such  a  scale  where  viewing  the  VFT  throught  the 
small  available  aperture  becomes  difficult.  To  overcome  this  limitation 
the  displacement  property  of  the  Fourier  transform  can  be  utilized. 

A  composite  transparency  containing  an  ordered  or  random  array  of 
reduced  replicas  of  the  transmittance  function  t(xQ,yo)  arranged  side 

by  side  as  illustrated  in  Fig.  4  is  prepared.  When  such  a  composite 
transparency  is  viewed  with  the  point  source,  the  VFT's  formed  by  the 
individual  elements  will  overlap  in  the  virtual  Fourier  plane.  The  VFT's 
are  identical  except  for  Linear  phase  dependence  on  x,y  which  depends 
in  each  VFT  on  the  central  position  of  each  element  in  the  composite 
transparency.  This  leads  to  a  desirable  noise  averaging  effect  and 
the  appearance  of  fine  checkered  texture  in  the  image  detail.  All 
this  leads  to  an  enhancement  of  the  quality  of  the  observed  power 
spectrum.  Both  coherent  and  noncoherent  superposition  of  the  over¬ 
lapping  VFT's  is  possible  using  this  scheme  by  varying  the  coherence 
area  of  the  wavefield  illuminating  the  composite  transparency.  When 
the  coherence  area  is  roughly  equal  to  the  size  of  the  individual 
elements  of  the  composite  transparency  noncoherent  superposition 
results,  while  a  coherence  area  equal  or  greater  than  the  size  of  the 
composite  transparency  would  yield  coherent  superposition. 


7 


Virtual  Fourier 
transform  plane 


Object  plane 
transparency 


Fig.  4.  A  composite  screen  consisting  of  an  ordered  array  of  identical 
Fourier'  transform  projection  holograms. 


THREE  DIMENSIONAL  DISPLAY 


The  VFT  concept  and  the  “weighted  Fourier  domain  projection  theorem" 
discussed  above  can  be  combined  in  an  attractive  scheme  for  the  recon¬ 
struction  and  display  of  a  3-D  image  from  a  series  of  weighted  projection 
holograms  corresponding  to  different  parallel  slices  through  the  object. 
The  scheme  is  based  on  viewing  a  series  of  weighted  projection 
holograms  sequentially  in  the  proper  order  of  the  occurance  of  their 
corresponding  slices  in  the  original  object  while  displacing  the  point 
source  axially  for  one  hologram  to  next  by  an  axial  increment  proportional 
to  the  spacings  between  adjacent  object  slices.  In  this  fashion  the 
reconstructed  virtual  images  of  the  various  slices  are  seen  in  depth 
at  different  VFT  planes  that  are  determined  by  the  positions  of  the 
axially  incremented  point  source.  Repeated  rapid  execution  of  this 
procedure  by  displacing  the  point  source  back  and  forth  leads  the 
observer  to  see  a  virtual  3-D  image  tomographically  in  parallel  slices 
or  sections  as  he  looks  through  the  series  of  projection  holograms 
passed  rapidly,  as  in  a  motion  picture  film,  infront  of  his  eyes. 

More  specifically  the  scheme  is  based  on  preparing  a  series  of 
N  weighted  Fourier  domain  projection  holograms  from  the  3-D  Fourier 
domain  data  F(w)  of  a  given  object  f(r)  as  described  in  the  preceeding 
sections.  Each  of  the  projection  holograms  would  correspond  to  a 
different  parallel  slice  through  the  object.  A  composite  transparency 
similar  to  that  shown  in  Fig.  4  is  formed  for  each  projection  hologram. 

In  fact  Fig.  4  is  an  example  of  a  computer  generated  composite  hologram 
containing  an  array  of  identical  weighted  projection  holograms  corres¬ 
ponding  to  one  slice  of  the  test  object  shown  in  Fig.  5.  The  test 
object  chosen  consisted  of  eight  point  scatterers  arranged  as  shown. 

The  3-D  Fourier  space  of  this  test  object  was  accessed  in  a  computer 
simulation  of  wavelength  diversity  imaging  as  described  in  a  companion 
paper  in  this  volume*. _  The  resulting  computer  generated  Fourier 
space  data  manifold  F(wl  was  used  to  compute  three  weighted  projection 
holograms  corresponding  to  the  three  planes  R^=]m,0,lm  of  Fig.  5 

containing  the  three  different  distributions  of  point  scatterers.  A 
composite  array  such  as  that  of  Fig.  4  was  formed  and  displayed  by  the 
computer  for  each  of  the  three  projection  holograms,  each  was  photo¬ 
graphed  yielding  a  set  of  three  projection  hologram  composite  trans¬ 
parencies.  Copies  of  these  were  then  mounted  on  a  rotating  wheel  as 
shown  in  Fig.  6  (a)  and  viewed  with  an  axially  scanned  point  source. 

Four  sets  of  transparency  copies  of  these  three  composite  projection 
holograms  were  mounted  in  the  order  1,2, 3, 2, 1,2  ...  on  the  periphery 
of  a  rotating  wheel  as  shown  in  Fig.  6  (a).  The  wheel  is  driven  by 
a  computer  controlled  stepper  motor.  The  axially  scanned  point  source 
was  produced  by  scanning  a  focused  laser  beam  back  and  forth  on  a 
length  of  fine  nylon  thread  with  the  aid  of  a  deflecting  mirror  mounted 
on  the  shaft  of  a  second  computer  controlled  stepper  motor  as  shown  in 
Fig.  6  (b).  The  laser  and  optical  bench  arrangement  for  forming  the 
scanned  focused  beam  appear  in  the  background  of  F.iq.  6  (a).  The 
computer  controlled  steppers  enable  precise  positioning  of  the  secondary 
point  source  on  the  scattering  thread  in  synchronism  with  the  hologram 

♦See  paper  entitled  "Holography,  Wavelength  Diversity  Inverse 
Scattering"  in  this  volume. 


Fig.  6. 


Quasi  real-time  three-dimensional  image  reconstruction  and 
tomographic  display  in  successive  slices  from  a  series  of 
projection  holograms  mounted  on  rotating  wheel  seen  in  fore 
front  of  (a);  Detail  of  laser  scanner  used  to  produce 
linearly  scanned  point  source  is  shown  in  (b). 


Photographs  of  three  slices  of  the  virtual  3-D  image  of  the  test  object  ot  ng.  d  oDta 
photographing  the  VFT's  formed  from  corresponding  Fourier  domain  projection  holograms. 
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being  viewed  so  that  the  VFTsare  formed  in  their  proper  planes.  A 
viewer  looking  at  the  axially  displaced  point  source  through  each 
transparency  mounted  on  the  wheel  as  it  passes  infront  of  his  eye 
will  see  a  3-D  virtual  image.  Photographs  of  the  three  virtual  images 
seen  by  an  observer  in  this  fashion  are  shown  in  Fig.  7.  An  opto- 
digital  scheme  for  rapid  real-time  implementation  of  the  procedure  realized 
above  is  shown  in  Fig.  8.  This  scheme,  presently  under  study,  utilizes 
a  rapid  recyclable  spatial  light  modulator  (SLM)  such  as  the  Itek  PROM 
in  order  to  form  VFT's  of  the  projection  holograms  displayed  by  the 
computer  in  real-time. 


CONCLUSIONS 

We  have  presented  the  basic  principles  of  tomographic  3-D  image 
display  based  on  Fourier  domain  projection  theorems.  One  possible 
implementation  of  the  principle  using  the  virtual  Fourier  transform 
and  a  series  Fourier  domain  projection  holograms  has  been  described. 

There  are  several  advantages  for  using  the  VFT  rather  than  the  real 
Fourier  transform  (RFT),  the  most  important  of  which  is  the  ease  with 
which  the  position  of  the  VFT  plane  can  be  moved  axially  by  simply 
moving  the  position  of  the  reconstruction  point  source.  The  VFT  approach 
was  adopted  in  the  present  study  because  it  is  much  easier  to  move  a 
point  source  rapidly  than  to  move  the  display  screen  needed  in  the  RFT 
approach.  Furthermore  focusing  in  the  VFT  approach  is  carried  out  by 
the  observer  while  in  the  RFT  approach  it  must  be  performed  by  the 
system.  Other  attractive  features  of  the  VFT  are: 

(a)  Simplicity  -  enables  direct  viewing  of  the  power  spectrum 
of  a  transparency  or  a  hologram  with  a  yariety  of  simple  point 
sources. 

(b)  The  scale  of  the  observed  VFT  can  be  easily  altered  by 
changing  the  distance  between  the  projection  hologram  transparency 
and  the  reconstruction  point  source. 

(c)  Lower  speckle  noise  and  therefore  higher  reconstructed  image 
quality  can  be  attained  by  using  nonlaser  point  sources  in  the 
reconstruction  such  as  LED  or  miniature  spectrally  filtered 
incandescent  lamps.  Further  reduction  in  speckle  noise  occurs 
when  an  array  of  the  projection  hologram  rather  than  a  single 
hologram  is  used  and  when  the  hologram  is  slightly  vibrated  or 
is  in  motion  because  of  a  noise  averaging  effect. 

(d)  Coherent  and  noncoherent  superposition  of  VFT's  is  possible 
by  altering  the  coherence  area  of  the  reconstruction  wavelfield. 

(e)  Because  of  the  Fourier  transform  nature  of  the  projection 
holograms  utilized,  the  resolution  requirements  from  the 
storage  medium  (photographic  film  or  the  CRT/SLM  system  of  Fig.  8) 
are  much  lower  than  would  be  needed  in  the  recording  of  a  Fresnel 
hologram  of  the  object  as  a  means  of  3-D  image  display.  The  3-D 
image  detail  contained  in  the  single  Fresnel  hologram  is  now 
distributed  over  a  series  of  lower  resolution  projection  holograms 
which  are  used  to  form  the  3-D  image  sequentially  in  time  in 
individual  slices. 


REFLECT 

MODE 

SLM 


digital  scheme  for  the  reconstruction  and  display  of  3-D  images  using  a  recyclable 
al  light  modulator  and  a  point  source  to  view  the  VFT  in  real-time. 
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(f)  Because  3-D  image  reconstruction  is  tomographic  (in  separate 
slices)  there  is  no  interference  between  the  wavefields  forming 
the  various  slices. 

(g)  Permits  other  forms  of  3-D  image  display  involving  spatial  or 
angular  multiplexing  in  a  fashion  similar  to  integral  holography. 
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This  thesis  investigates  the  development  of  a  broadband 
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phase  shift  on  holographic  frequency  diversity  imaging  is 
considered  and  techniques  for  the  removal  of  this  phase 
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holograms  are  compared  to  those  measured  experimentally. 
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I  INTRODUCTION 


Frequency  diversity  imaging  has  been  under  study  at  the 
Electro-Optics  and  Microwave  Optics  laboratory  of  the  Moore 
School  Graduate  Research  Center.  [I  ],[2],[3],[<i]  This 
study  has  established  the  theoretic  feasibility  of  imaging 
objects  by  means  of  their  multiaspect  frequency  response  . 
For  the  purpose  of  experimentally  studying  frequency 
diversity  imaging,  an  experimental  measurement  system  has 
been  assembled  and  installed  in  the  Moore  School  anechoic 
chamber.  In  this  thesis  we  will  describe  the  automation  of 
this  measurement  system  and  characterize  its  peroformance  in 
the  measurement  of  complex  field  amplitudes  of  scattered 
fields.  A  system  block  diagram  fig. 1.1,  shows  the  major 
system  components.  The  central  element  is  the  DEC  MINC  LSI 
11/2  minicomputer.  This  computer  performs  several  important 
functions.  The  MINC  controls  laboratory  instrumentation  via 
the  IEEE-488  bus  protocol  standard.  This  allows  the 
Hewlett-Packard  8620C  microwave  sweeper  to  be  precisely 
tuned  to  any  frequency  in  the  2.0  to  18.0  GHz  range 
(fig. 1.2).  The  computer  collects  data  from  the  HP  8410B 
network  analyzer  through  the  four  analog  input  channels 
available.  These  analog  values  are  proportional  to  the 
amplitude  in  db  and  phase  in  the  range  (-TT)  to  (+ff)  radians 
relative  to  the  reference  signal  supplied  to  the  network 
analyzer.  The  computer  stores  Experimental  data  on  floppy 
discs.  The  available  storage  capacity  is  large,  over 


500,000  measurement  pairs  of  complex  field  amplitude  and 

I 


Fig.  1.1  System  block  diagram. 


if. 


Fig.  1.2  HP  8620C  microwave  sweeper  and  HP  1410B  network 
analyzer . 
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phase  may  be  stored  on  a  single  disc.  This  data  can  be 
accessed  for  both  processing  and  display  on  a  Tektronix  606A 
high  resolution  CRT  monitor.  The  disc  system  allows  the 
MINC  to  operate  under  a  sophisticated  software  system,  DEC 
RT-11  V3.0B  permitting  programing  in  MINC  BASIC, FORTRAN  IV, 
and  MACRO  languages. 

The  processing  capability  of  the  system  allows  the 
removal  of  system  response  errors  due  to  anechoic  chamber 
clutter , antenna  cross-coupling  and  receiver  channel 
characteristics.  The  data  may  be  processed  for  target  range 
calculations  and  the  removal  of  the  phase  shift  due  to  the 
target  range.  In  addition  the  collected  data  may  be 
filtered  to  improve  imaging  and  finally  displayed  on  a  high 
resolution  Tektronix  606A  X-Y  CRT  monitor  using  the  MINC 
system  D/A  converter  module. 

The  second  part  of  this  thesis  will  describe  the 
simulation  and  actual  operation  of  a  holographic  radar 
system  verifying  the  theory  of  frequency  diversity  imaging. 
This  is  done  using  the  system  described  in  the  first 
section.  The  scattering  of  various  targets  is  derived  and 
holograms  using  these  results  are  generated  for  comparison 
with  experimental  data.  Finally  a  system  for  the 
experimental  measurement  of  the  scattering  for  these  targets 
is  outlined  and  the  results  from  this  system  compared  to 
theory. 
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II  SYSTEM  OPERATION  AND  ERROR  REMOVAL 

This  chapter  will  cover  the  operation  of  the  microwave 
backscatter  data  aquisition  system.  This  will  include 
actual  interfacing  information  and  an  analysis  of  the  types 
of  errors  encountered  when  making  microwave  measurements. 
The  error  correction  techniques  developed  are  later  used  in 
the  experimental  verification  of  the  frequency  diversity 
imaging  theory. [ H  ] 

2.1  Microwave  sweeper  operation 

The  first  task  in  the  development  of  an  automated  data 
acquisition  system  is  the  implementation  of  a  data 

communications  link  between  the  intelligent  controller  and 
instrumentation.  The  IEEE-488  bus  protocol  is  utilized  in 
this  application  for  the  transfer  of  data  to  the 

Hewlett-Packard  8620C  microwave  sweeper  from  the  MINC  LSI 
11/2  computer.  This  bus  is  a  high  speed  8  bit  wide 

bidirectional  data  path  with  5  additional  lines  dedicated  to 
control.  Data  is  transferred  in  ASCII  format  over  the  bus. 
For  example  the  number  1  is  transmitted  as  the  ASCI  code  for 
the  character  'l'.  Certain  sequences  of  characters  make  the 
sweeper  perform  different  functions  or  enter  different  modes 
of  operation  via  its  IEEE-488  interface. 

In  order  to  set  the  frequency  of  the  sweeper  a  number 
must  be  sent  to  the  IEEE-488  interface  in  the  range  1-10000. 
Each  sweeper  frequency  band  has  been  split  into  10000 
frequency  points.  The  frequency  of  operation  is  controlled 


by  an  internal  analog  voltage  that  varies  between  0.0  ana 

10.0  volts.  A  D/A  converter  on  the  8620C  IEEE-488  interface 

changes  the  data  transmitted  from  the  MINC  into  the 

frequency  controlling  voltage.  The  correspondence  between 

voltage  and  output  frequency  is  essentially  linear.  In 

order  to  obtain  the  interpolating  function  for  frequency 

versus  control  voltage;  a  microwave  frequency  counter  was 

used  to  measure  the  voltage-frequency  characteristic 

function.  Using  the  MINC-BASIC  program  CALAB.BAS  a  least 

squares  fit  for  both  linear  and  quadratic  functions  was  made 

on  the  frequency  vs  voltage  data.  This  type  of  program  is 

used  for  determining  the  best  polynomial  fit  to  the 

f requency-voltage  characteristic  of  the  sweeper.  Sample 

output  and  program  listings  are  in  appendix  I  along  with  an 

explanation  of  program  operation.  The  results  from  this 

work  indicate  that  the  quadratic  fit  was  statistically 

superior  for  all  bands  on  the  microwave  sweeper.  The 

FORTRAN  subroutine  SWEEP  was  written  which  utilizes  the 

quadratic  interpolation  polynomial  for  each  of  the  four 

bands  of  the  8620C.  This  subroutine  automatically 

calculates  the  control  voltage  and  band  to  generate  any 

frequency  in  the  2-18  GHz  range  and  transmits  the 

appropriate  commands  over  the  IEEE-488  bus.  The  variance  of 

the  frequency  setting  using  the  quadratic  fit  for  the  three 

bands  are  as  follows:  .6  MHz  in  the  2.0-6. 3  GHz  band;  1.2 

MHz  in  the  6.3-12.0  GHZ  band;  and  1.6  MHZ  in  the  12.0-18.0 

GHz  band.  For  higher  accuracy  the  sweeper  may  be  phase 
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locked  to  the  reference  in  a  locking  frequency  counter 
yielding  very  high  accuracy  as  precise  as  the  frequency 
reference  itself.  The  EIP  371  locking  counter  may  be  used 
in  this  application  to  lock  the  HP  8620C  sweep  oscillator  to 
the  correct  frequency  once  it  is  within  20  MHz  of  the 
desired  frequency.  The  auxiliary  output  of  the  HP  8620C 
supplies  a  sample  of  the  signal  generated  by  the  fundamental 
2. 0-6. 3  GHz  oscillator  module  within  the  sweeper  to  the 
locking  counter.  Sweeper  output  on  higher  bands  is  this 
fundamental  multiplied  by  a  factor  of  2  or  3  for  the 
6.3-12.0  and  12.0-18.00  GHz  bands  respectively.  For  example 
the  locking  frequency  for  a  10.0  GHz  sweeper  output  on  band 
2  would  be  5.0  GHZ.  In  operation,  subroutine  SWEEP  will  set 
the  frequency  of  the  sweeper  within  2.0  MHz  of  the  desired 
frequency  and  subroutine  SLOCK  will  be  called  to  calculate 
the  locking  frequency;  lock  the  HP  8620C  and  return  to  the 
main  calling  routine  when  lock  will  have  occurred.  Lock 
time  varies  from  .1  to  3  seconds  and  resolution  is  100  KHZ. 
These  subroutines  are  called  whenever  the  sweeper  must  be 
set  to  a  particular  frequency  or  the  sweeper  must  be  placed 
in  or  be  released  from  computer  control. 

2 . 2 _  Ne twork  analyzer-compute r  interface 

The  Hewlett-Packard  8410B  network  analyzer  is  the  focus 
of  the  measurement  capabilities  of  the  microwave  measurement 
system.  It  can  make  vector (amplitude  and  phase) 
measurements  in  the  (.l-18.0)GHz  range.  The  range  of 
amplitude  measurement  is  80^  db  and  phase  may  be  measured 


modulo  (2TT)  .  The  system  reference  signal  is  fed  from  a  20db 
directional  coupler  to  the  HP8411A  harmonic  converter 
sampling  head  of  the  network  analyzer.  This  reference  is 
compared  to  the  backscatter  from  the  illuminated  target; 
both  in  amplitude  and  phase.  The  reference  signal  amplitude 
is  kept  constant  by  leveling  the  sweeper  with  a  feedback 
signal  derived  from  its  amplified  output  by  means  of  a 
crystal  detector.  This  allows  the  sweeper-TWT  (Traveling 
Wave  Tube)  system  to  yield  nearly  constant  output  in  the 
(2.0-16.5)  GHz  range;  see  fig. 2.1.  TWT  power  output  is  on 
the  order  of  1  watt  over  these  frequencies. 

The  complex  field  amplitude  measurements  are  available 
as  analog  voltages  from  the  back  panel  of  the  8410B.  The 
outputs  are  proportional  to  amplitude  and  phase:  25  MV/db 
and  10  MV/DEGREE  These  values  are  digitized  by  the  MINC 
using  its  built  in  analog  to  digital  conversion  channels. 
The  MINC  A/D  converters  digitize  voltages  lying  in  the  range 
of  -5.12  to  +5.12  volts  to  the  range  of  0-4096  yielding  12 
bit  resolution.  If  the  signal  is  corrupted  by  noise;  the 
user  has  the  option  of  employing  signal  averaging  to  cancel 
the  effects  of  noise  uncorrelated  to  the  received  signal. 

A  difficulty  encountered  when  measuring  phase  angle 

modulo  ( 2tt)  occurs  when  the  phase  is  close  to  (+ir)  or  (-rr)  . 

At  this  point  a  small  change  in  the  signal  phase  may  cause 

the  phase  to  flip  between  these  two  equivalent  extremes 

rapidly.  If  a  data  sample  is  taken  close  to  (+tr)  or  (— rr)  it 

may  be  in  transition  between  them  and  therefore  incorrect. 
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Since  such  points  occur  infrequently  in  a  typical 
measurement  they  may  be  ignored.  ; their  presence  does  not 
seriously  hinder  any  holographic  imaging  due  to  the  inherent 
redundancy  and  therefore  noise  immunity  of  the  holographic 
reconstruction  process. [5]  If  it  is  desired  that  these 
points  be  identified  and  removed  it  is  necessary  to  estimate 
the  mean  and  variance  of  the  samples  taken  at  each  frequency 
point.  At  frequencies  where  the  phase  is  flipping  between 
(+tt)  and  (-TT  )  ,  the  variance  will  be  much  larger  than  at 
other  frequencies.  The  mean  of  the  samples  when  this  is 
occurring  will  be  near  zero.  Hence  to  resolve  the  ambiguity 
problem  it  is  necessary  to  decide  if  the  variance  exceeds  a 
predetermined  threshold  when  the  mean  in  the  neighborhood  of 
zero.  Two  FORTRAN  subroutines  PHAMP  and  PHAMP2  which 
implement  these  algorithms  are  listed  in  appendix  I. 
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2.3 


Implementation  of  the  data  aquisition  system 


In  an  experimental  environment  it  is  important  to  be 
aware  of  the  various  types  of  errors  inherent  in  the 
equipment  and  the  experimental  procedure  adopted. 
Conditions  and  equipment  always  vary  from  theoretical 
ideals.  A  clear  understanding  of  the  error  removal  process 
leads  to  the  development  of  practical  implementations  of 
theoretical  concepts  and  enhancement  of  measurement  accuracy 
unattainable  otherwise. 

Errors  in  complex  field  amplitude  measurement  may  be 

caused  by  several  factors.  These  may  be  grouped  into  two 

categories.  Errors  caused  by  the  instruments  themselves 

fall  into  the  first  class.  Such  factors  as  measurement 

variations  caused  by  electronic  noise,  , inaccurate  A/D  and 

logarithmic  conversions,  and  inaccuracy  and  instability  in 

the  microwave  source  make  up  this  category.  The  second 

group  of  errors  is  caused  by  the  test  set, antennas  ,  cables, 

connectors,  amplifier,  and  room  clutter.  All  these  factors 

interact  with  each  other  in  the  microwave  region  and  are  the 

significant  cause  of  error  in  microwave  measurements. 

Little  can  be  done  about  the  first  class  of  errors  since 

they  are  inherent  in  the  characteristics  of  the  equipment 

used.  The  second  group  of  errors  can  be  removed  through  the 

use  of  automated  measurement  of  system  parameters  in  the 

frequency  range  of  interest.  These  errors  can  be  removed 

form  any  measurements  of  scattering  objects  by  digital 

processing  and  the  results  stored  for  later  recall.  This 
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essentially  provides  an  automated  and  improved  version  of 
the  conventional  two  antenna  radar  cross-section  measurement 
technique  [6]  ;  in  which  a  microwave  bridge  is  balanced  in 
the  absence  of  the  target  and  the  degree  of  imbalance  is 
measured  when  the  target  is  introduced  into  the  microwave 
field. 

Let  us  look  at  the  first  class  of  errors  more  closely 
since  these  errors  will  set  the  ultimate  performance  limits 
on  the  system.  The  characteristics  of  the  signal  source  are 
important  in  this  regard.  In  this  case  the  signal  source  is 
a  Hewlett-Packard  8620C  microwave  sweeper.  Since  the 
sweeper  is  not  phase  locked  frequency  and  stability  'problems 
exist.  The  carrier  also  has  significant  FM  noise  which 
appears  as  phase  noise  in  the  scattered  signal.  The  phase 
shift  of  the  scattered  scattered  signal  as  a  function  of 
frequency  for  small  frequency  variations  is  given  by: 

AG  =  A  k-R  u.i) 

where  k=(2TT/*)  and  R  is  the  path  length.  For  R  greater  than 
a  few  tens  of  meters  ,the  FM  noise  on  the  signal  source 
causes  measurable  variations  in  the  phase  of  the  scattered 
signal.  The  stability  of  a  synthesizer  is  required  for  the 
implementation  of  a  holographic  radar  system  when  target 
ranges  are  in  terms  of  kilometers. 

Another  limit  on  the  ultimate  accuracy  of  the  system  is 

the  resolution  of  the  A/D  conversions  and  the  accuracy  of 

the  network  analyzer.  Given  the  2.2MV  resolution  of  the 

M1NC  A/D  converter  and  the  network  analyzer  analog  output  of 

/£> 


50  MV/db,  the  system  can  resolve  .048  db  steps.  This 

resolution  limit  restricts  the  minimum  signal  to  system 

error  ratio.  If  the  error  signal  consisting  of  clutter, 

antenna  coupling,  system  directivity  and  noise  exceeds  the 

scattered  target  signal  the  resolution  of  the  target  signal 

suffers.  For  example,  if  the  system  error  signal  and  target 

scattered  signal  are  of  equal  intensity,  then  a  1  db  change 

in  the  scattered  signal  causes  a  .53  db  change  in  the  total 

received  signal  vector.  When  the  system  error  is  13  db 

above  the  scattered  signal  it  becomes  impossible  to  resolve 

a  1  db  change  in  the  target  signal  given  the  resolution 

capabilities  of  the  system.  This  difficulty  is  further 

compounded  by  errors  in  the  network  analyzer.  These  too  are 

amplified  when  the  clutter  exceeds  the  target  scatter 

signal.  In  fig. 2.1  is  a  plot  of  the  minimum  system 

resolution  in  order  to  detect  a  1  db  change  in  the  target 

return  signal  versus  the  system  error  signal  to  scattered 

signal  level  in  db.  In  fig. 2. 2  is  plotted  the  target  signal 

resolution  versus  the  noise  to  signal  level  in  db.  When  the 

system  error  is  20  db  below  the  target  signal  then  the 

resolution  of  the  target  scatter  signal  is  very  close  to  the 

ultimate  resolution,  .048db.  Clutter  is  the  component  of 

the  received  signal  not  scattered  by  the  target  but  that 

signal  that  is  the  result  of  coupling  between  antennas  and 

signal  scattered  by  the  anechoic  chamber  walls.  As  the 

clutter/signal  ratio  increases  the  target  scatter  signal 

resolution  decreases  exponentially.  Clutter  may  be  reduced 
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Fig.  2.1  The  required  resolution  in  db  for  a  data 
acquisition  system  to  detect  a  1  db  change  in  the  desired 
signal  vs.  signal/error  -signal  ratio. 
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Fig.  2.2  Minimum  change  in  desired  signal  level 
versus  signal/error  signal  level  given  .048  db 
acquisiton  resolution. 


detectable 
system  data 


by  improving  the  isolation  between  the  transmitter  and 
receiving  antennas  with  the  introduction  of  absorbing  foam 
panels  such  as  Emerson  and  Cummings"  Ecosorb  panels  between 
the  antennas. 

2.4  System  error  correction 

Turning  next  to  the  second  class  of  errors;  those 
directly  measurable  and  therefore  removable;  define  the 
following  quantities  which  are  functions  of  frequency  ; 

1(f)  —  Isolation  of  reference  to  test  channels 
T(f)  —  Transfer  characteristic  of  system 
A (f )  —  Attenuator  characteristic 
Al(f) —  Antenna  system  characteristic 
S(f)  —  Corrected  backscatter  for  target 
Cl(f) —  Uncorrected  antenna  clutter  and  coupling 
C2(f) —  Corrected  antenna  clutter  and  cross-coupling 
(uncorrected  for  antenna  system  response) 
C(f)  —  Corrected  antenna  clutter  and  coupling 
Rl(f) —  uncorrected  reference  target  backscatter 
R2(f) —  Corrected  reference  target  data 

(uncorrected  for  antenna  system  response) 
R(f)  —  Corrected  reference  target  data 
Several  possible  techniques  exist  for  the  removal  of  system 
errors.  The  particular  technique  is  dependent  on  the 
relative  signal  levels  involved,  the  accuracy  desired,  ease 
of  implementation,  and  computational  speed.  The  first 
technique  described  here  is  similar  to  that  used  by  Weir  et 
al  .(7] 
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The  first  step  in  the  correction  procedure  measures  the 
transfer  function  of  an  attenuator  A(f)  as  a  function  of 
frequency.  The  equipment  setup  for  this  procedure  is  shown 
in  fig.  2.3.  The  two  ports  of  the  reflection-transmission 
unit  connected  through  a  precision  HP  11605A  flexible 
coaxial  arm.  The  MINC  then  steps  the  sweeper  to  a  number  of 
frequency  points  and  stores  the  system  response  (log 
amplitude  and  phase  vs.  frequency)  in  memory.  When  this 
completed  the  attenuator  is  placed  in  series  with  the  arm 
and  another  set  of  measurements  is  made  at  the  same 
frequency  points  as  before.  the  system  response 
characteristic  is  subtracted  from  the  combined  attenuator 
plus  system  response  measurement  made  on  the  second  sweep. 
An  example  of  this  procedure  is  shown  in  fig  2.4.  Computer 
subroutine  PAD  performs  this  operation.  It  is  listed  in 
appendix  I  along  with  all  other  computer  program  listings 
and  output  pertaining  to  system  response  measurement  and 
removal. 

The  next  step  in  the  calibration  process  is  measurement 

of  the  reference  to  test  channel  isolation  I (f )  .  This 

characteristic  is  dependent  on  the  directivity  of  the 

network  analyzer  harmonic  converter , and  the  reflection 

transmission  unit.  For  the  Hewlett-Packard  8411A  harmonic 

converter  the  isolation  is  greater  that  50  db.  In  this 

measurement  the  ports  of  the  reflection-transmission  unit 

are  terminated  in  the  cables  used  for  the  later  target 

scatter  measurement.  These  coaxial  cables  are  terminated  in 

IH 


Fig.  2.4  30  db  attenuator  characteristic  for  8.0-8. 5  GHz. 
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50  ohm  resistive  loads  as  shown  in  fig. 2.5.  The  results  of 
a  typical  run  are  shown  in  fig. 2.6.  Computer  subroutine  1ST 
automates  this  stage  in  the  correction  procedure.  As  can  be 
seen  the  coupled  signal  is  well  in  the  noise  of  the  system 
and  would  not  affect  later  scatter  data.  If  the  isolation 
effect  is  ignored  later  calculation  would  be  simplified 
greatly;  but  is  included  here  to  be  consistent  with  the 
procedure  outlined  in  the  literature. [7] 

The  next  stage  in  generating  the  data  for  correction  of 
scatter  data  is  measurement  of  the  system  transfer  function 
T{f).  This  consists  of  the  characteristics  of  the  traveling 
wave  tube  amplifier , system  cables  and  connectors.  This  is 
done  by  connecting  the  cable  from  the  transmitting  antenna 
to  the  receiving  antenna  cable  in  fig. 2.6  and  placing  the  30 
db  attenuator  characterized  previously  in  the  line  to  avoid 
damage  to  the  harmonic  converter.  The  raw  measurement  MT (f ) 
is  a  combination  of  several  factors: 

MT  ( f ) =T ( f ) * A ( f ) +1 ( f )  (2.2) 

Solving  for  T{f): 

T(f)*(MT(f)-I  (f)  )/A(f)  (2.2a) 

The  equipment  setup  for  this  procedure  is  shown  in  fig  2.7 
and  typical  uncorrected  and  corrected  transfer  function  data 
in  fig. 2.8  and  2.9. 

Measurement  of  the  antenna  cross-coupling  and  room 

clutter  is  the  next  step  in  the  correction  process.  In  this 

procedure  shown  in  fig. 2. 7b,  the  target  is  removed  from  the 

anechoic  chamber  and  the  antennas  pointed  to  the  target 
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Pig.  2.5  System  configuration  for  measuring  isolation 
between  reference  and  unknown  ports. 
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Pig.  2.6  Isolation  between  ports  of  reflection-transmission 
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Fig.  2.7  a)  System  configuration  for  transfer  characteristic 
measurement.  b)  Connection  to  antennas  for  scattering  and 
clutter  measurement.  c)  Transmitting  antenna  with  spiral 
AEL  antenna  and  parabolic  dish.  d)  Receiving  dual 
polarization  horn  ;  EMI  A6100  with  Norsal  90  hybrid. 


Fig.  2.7  (contd.)  a)  system  configuration  for  transfer 
characteristic  measurement.  b)  Connection  to  antennas  for 
scattering  and  clutter  measurement  c)  Transmitting  antenna 
with  spiral  AEL  antenna  and  parabolic  dish.  d)  Receiving 
dual  polar iztion  horn;  EMI  A6100  with  Nor sal  90  hybrid. 
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Fig.  2.8  Uncorrected  transfer  characteristic 
8. 0-8. 5  GHz. 


Fig.  2.9  Transfer  characteristic  corrected  for 
response;  8.0-8. 5  GHz. 
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location.  A  high  gain  parabolic  dish  antenna  is  used  for 
illumination  of  the  target  since  the  narrow  beam  pattern  of 
the  antenna  places  most  of  the  radiated  power  on  the  target 
area.  A  smaller  dual  polarization  horn  is  used  for 
receiving  in  order  to  sample  a  small  area  of  the  scattered 
field.  The  uncorrected  clutter  Cl(f)  is  given  by: 

Cl(f)=C2(f)*T(f)+I  (f)  (i  s) 

Solving  for  the  corrected  clutter  and  coupling: 

C2(f)  =  (Cl  ( f )  — I  (f)  )  /T  ( f )  (2.3  cS) 

Subroutine  ANTEN  does  the  system  clutter  and  coupling 
removal.  Examples  of  the  uncorrected  and  corrected  clutter 
Cl(f)  and  C2(f)  are  shown  in  figs. 2. 10  and  2.11.  The 
corrected  clutter  represents  the  actual  signal  reflected 
from  the  anechoic  chamber  walls  and  that  signal  coupled 
between  the  antennas  with  the  system  response  removed. 

These  subroutines:  PAD , 1ST, TRANS ,  and  ANTEN,  were 
combined  into  a  program  SYSRES.  Data  from  each  of  these 
subroutines  may  be  stored  on  disc  for  later  recall  or 
display.  Theoretically  if  the  system  is  not  disturbed  then 
the  system  response  will  remain  constant.  Then  only  the 
target  data  need  be  recorded  in  any  run  for  a  new  corrected 
backscatter  measurement. 

The  transfer  function  of  the  system  and  the  range 

clutter-  antenna  cross  coupling  data  are  utilized  the  the 

next  step  of  the  error  correction  process.  A  reference 

object  of  known  constant  cross  section  is  measured  and  the 

result  stored.  This  data  includes  all  the  errors  previously 
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Fig.  2.11  System  clutter  corrected  for  transfer  function; 
8. 0-8. 5  GHz. 
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described  but  also  takes  into  account  the  antenna  system 


variations  as  a  function  of  frequency: 

Rl<f  )=C2(f  )*T(f  )+I (f )+R2(f)*T(f ) 


(Z.M) 


Here  C2(f)  and  R2(f)  contain  the  antenna  system  response 
multiplying  the  actual  values  of  corrected  reference  target 
and  system  clutter  data. 

R2(f )=R(f)*Al(f)  (2.S) 

C2  (f )  =C  (f )  *A1  (f )  (2.fc) 

This  response  Al(f)  takes  into  account  the  varying  amount  of 
power  received  and  transmitted  as  a  function  of  frequency  in 
the  antenna  system.  If  the  reference  target  is  chosen  to 
have  a  constant  cross  section  and  linear  phase  over  the 
frequency  range  of  interest  then  R(f)  is  of  constant 


amplitude  and  linear  phase.  Solving  for  R2(f): 
R2 (f)=(Rl (f)-I (f) -Cl (f ) )/T (f ) 


C2.7) 


This  leaves  R2(f)  proportional  to  Al(f)  shifted  by  a  linear 

phase  corresponding  .to  the  reference  target  range. 

When  the  actual  target  is  measured;  it  is  corrected  for 

system  errors  as  was  the  reference  target  data  and  this 

result  is  divided  the  corrected  reference  target  data  to 

yield  the  final  target  scatter  data. 

S2(f)  =  (Sl(f)-I  (f)-Cl(f) )  /T  (f )  (2.8) 

S  (f  )=S2  (f  )/R2  (f )  (Z  9) 

This  technique  may  be  simplified  considerably  in  the 

laboratory  environment  given  the  signal  to  noise  ratio  is 

greater  than  10  db  for  the  scattered  signals.  In  this  case 

only  multiplicative  errors  remain  and  the  additive  errors 

Zh 


are  masked  by  the  high  target  scatter  signal  amplitude. 
Hence: 


1(f)  ,Cl(f)->0 
R1 (f ) =R2 (f ) *T(f) 


(2.10) 


and  therefore: 

S  (f )  =S1  (f  )/Rl  (f )  (2.11) 

Weir  and  his  group  have  reported  that  this  technique  has 
reduced  equivalent  range  clutter  to  -45  db  below  1  sq. 
meter . 

There  remains  one  source  of  error  in  the  scattered 
signal  measurement  that  cannot  be  removed  by  calculation. 
This  error  comes  from  multipath  scattering  from  the  object. 
The  target  scatters  power  in  all  directions.  Some  of  this 
signal  may  be  reflected  off  the  walls  or  floor  of  the 
anechoic  chamber.  The  signals  reflected  off  the  walls  is 
over  48  db  down  from  the  incident  wave  amplitude  in  the 
6.0-12.0  GHz  range.  However  if  the  target  is  small  in  cross 
section;  on  the  order  of  100  sq.  cm.;  then  it  is  possible 
for  the  walls  (on  the  order  of  100  000  sq.  cm.  )  to 
contribute  a  significant  component  to  the  received  signal. 

In  order  to  analyze  the  effect  of  the  multipath 
scattering,  let  the  directly  received  signal  be  written: 


SB(t)  *  A  eos(u»t  ♦  V) 


(2.11) 


and  the  indirectly  received  signal: 


S;(t)  *  B  c-©s(ujt  «y*e) 

Then  the  total  received  signal  is  given  by: 
S,(t)  *  («»♦  2A8<oS«»*B*)  coi(wt  ♦  y* 


U.I3) 


(MH) 
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Since  theta  is  a  function  of  the  indirect  path  length 
differences  and  frequency,  it  will  lead  to  a  periodic 
variation  in  the  amplitude  of  the  scattered  signal.  As  an 
example;  if  the  paths  differ  in  length  by  1  meter  then  the 
amplitude  oscillations  will  occur  every  300  MHz  given  all 
other  factors  remain  constant. 

A  series  of  programs  was  written  to  test  these  various 
techniques  of  error  removal.  They  differ  in  the  only  in  the 
error  removal  technique  employed;  not  in  file  storage  or 
display  formats  nor  in  range  calculation  and  removal  to  be 
described.  These  programs  used  together: 

1) Measure  and  correct  the  reference  target  data  with 
files  of  transfer  function  and  clutter  data 
generated  by  SYSRES. 

2) Measure  and  correct  the  -target  data  and  finally 
take  the  corrected  reference  target  data  and  remove 
the  antenna  system  response  from  the  target  data. 

3)  Alternately  for  high  SNR;  calculate  the  correct 
scattered  target  signal  directly  using  the  reference 
target  data. 

4)  Calculate  and  remove  phase  shift  due  to  range  from 
the  target  signal  after  error  correction 

These  programs  are  briefly  described  and  differ  in  the 
mentioned  categories: 


SPHERE-Reads  transfer  function  and  clutter  data 
files  generated  by  SYSRES;  takes  the  reference  or 
object  data  and  corrects  for  errors  in  the 


equipment. 

SPHER2-Measures  clutter  and  reference  target  data. 
It  then  subtracts  clutter  and  divides  the  target 
data  by  the  reference  target  transfer  function. 
SPHER3-Measures  reference  target  and  object  signals 
ignoring  clutter,  and  divides  the  object  data  by  the 
reference  target  complex  field  amplitude  data. 

In  order  to  implement  the  complete  error  correction  process 
program  SPHERE  would  be  run  twice;  once  for  the  reference 
target  and  then  again  for  the  test  object.  This  data  would 
be  stored  for  later  retrieval.  Program  SPHER3  would  read 
these  files  and  process  then  such  that  the  test  object 
complex  field  data  would  be  divided  by  the  reference  target 
data.  For  high  SNR  cases;  program  SPHER3  alone  would  be 
runs  first  measuring  the  reference  target  and  then  the  test 
target  and  finally  dividing  them  yielding  the  final  result. 
Listings  and  a  more  complete  description  of  program 
operation  is  given  in  the  program  appendix  I. 
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Ill  RANGE  PHASE  SHIFT  ANALYSIS 


This  chapter  contains  an  analysis  of  the  effects  on  the 
phase  shift  due  to  target  range  on  coherent  imaging  of  a 
target.  It  includes  a  discussion  of  some  of  the  available 
techniques  for  the  removal  of  this  phase  shift  and  the 
required  perfor  mance  of  these  systems  based  on  bandwidth 
and  signal  to  noise  ratio. 


3.1  The  effects  of  range  phase  on  imaging 

As  previously  described  by  Farhat  and  Chan  ; HI ;  the 
scattered  field  of  the  scattering  target  is  given  by: 


j"U  e-0 


Where  Rt  and  R,  are  the  vectors  from  the  receiving  and 
transmitting  antennas  to  the  target  respectively  and  k  is 
the  wave  vector  of  the  illuminating  wave.  In  this  case  the 
integration  is  over  the  extent  of  the  object.  The  integral 
term  is  independent  of  the  range  of  the  target.  While  the 
term  preceeding  the  integral  is  target  independent  and 
contains  range  information.  The  argument  of  the  complex 
exponential  is  a  linear  phase  function  of  frequency. 

In  the  single  receiver=transmitter  pair  arrangement; 

shown  in  fig. 3.1,  the  spatial  frequency  domain  (p  space) 

data  is  collected  in  a  plane  perpendicular  to  the  axis  of 

rotation.  In  this  case  the  multiplication  of  the  range 

phase  factor  leads  to  the  convolution  of  the  transforms  of 

the  the  range  and  target  functions  in  the  spatial  domain. 

The  real  part  of  the  range  frequency  domain  function  is  a 
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radially  symmetric  sinusoidal  function: 

fife  *  cos  -h  { Rt  <Qf)  (J.2) 

This  transforms  to  a  circular  ring  of  radius  (Rr+  iy  .  This 
indicates  that  each  point  of  the  reconstructed  image  is 
convolved  with  this  circular  ring  pattern;  seriously 
degrading  the  target  image.  In  fact  any  error  in  the 
removal  of  phase  will  distort  the  image  in  this  manner.  The 
effect  of  removing  the  range  phase  factor  is  equivalent  to 
the  focusing  of  the  system  on  the  target. 


There 
factor : 


are  several  other 


e 


reasons  for 

•ikto'*') 


the 


removal  of  the 
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from  the  received  data.  If  the  data  is  discrete  then  the 
considerations  of  aliasing  and  sufficient  data  sampling  rate 
are  introduced.  When  the  sampling  rate  in  the  frequency 
domain  is  (  f)  then  the  maximum  target  range  before  aliasing 
will  occur  is  (c/4*  f) .  However  if  the  phase  factor  is 
removed  and  the  target  is  of  smaller  dimension  L  than  the 
range,  then  the  sampling  rate  in  the  frequency  domain  need 
only  be  sufficient  to  prevent  aliasing  over  the  dimensions 
of  the  object  : 

(3V 


This  allows  the  entire  resolution  capability  of  the  system 

to  be  placed  on  the  target  itself  greatly  reducing  the  data 

2? 


m 


volume . 


Another  reason  for  the  removal  of  this  phase  factor 
term  can  be  seen  in  systems  involving  multiple 
receiver-transmitter  pairs.  For  each  pair  (|CT)  and  (Rt )  is 
different  and  in  order  for  the  data  to  be  coherent  all  phase 
centers  must  be  equal  for  an  image  to  be  formed. 

Several  methods  have  been  suggested  for  the  removal  of 
the  range  phase  factor.  [4]  In  all  cases  it  is  required 
that  the  range  removal  technique  be  accurate  to  within  (a/5) 
for  there  not  to  be  serious  image  degradation.  It  is 
important  to  investigate  the  constraints  on  the  ranging 
system  parameters  necessary  to  attain  the  required  accuracy. 

From  the  scaling  theorem  in  Fourier  analysis;  a  signal 
cannot  be  both  of  narrow  bandwidth  and  short  duration.  [8] 

*F(£)  (3.5) 

We  define  the  duration  and  the  bandwidth  of  the  signal  in 
the  following  manner: [  ] 

/At)*  *  J  /°t*|i(t)J*Jt  (a)  dua  ft,)  (3.fc«.,b) 

where: 

E = .  h£  IM  a  40 

is  the  signal  energy.  Then  if: 

Jr*  VT"  /ft)  (3.U) 

\t\*eo 

When  applyed  to  signals  scattered  by  a  target  this  may  be 

translated  to  range  uncertainty: 

AR  Li 
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(3.7) 


This  represents  a  limit  when  the  durations  of  the  signal 
pairs  are  as  previously  defined.  However  it  may  be  possible 
to  improve  the  range  resolution  given  that  the  signal  to 
noise  ratio  is  greater  than  10  db. 

Consider  a  sinusoid  in  narrow  band  gaussian  noise: 

r (t\  z  r  to*  ui  t  •*  s  ~  '''s  It}  ai*.  ( 3.f ) 

nW  *  ntW  - 'ftjU)  Suv  ui  "t 

where  n(t)  is  the  noise  signal  and  ns(t)  and  nc(t)  are  the 
quadrature  components  of  the  noise  signal.  The  amplitude  and 
phase  may  be  expressed  as: £9} 

l^l  :  ^  (3.9a.) 

J  (3  »t) 

Since  the  noise  is  gaussian  nc(t)  and  ns(t)  are  gaussian. 
Assuming  that  the  SNR  is  high  the  phase  may  be  approximated 


Also  since  the  noise  is  white  n  %  (t)  may  be  related  to  the 
bandwidth  of  the  receiving  system. 

<=*^B  (3.U) 

where  (•;*)  is  the  variance  of  the  n s(t)  quadrature  phase 
component.  The  probability  density  of  the  phase  may  be 
written  as: 

A  i 

c  *Nf®.  fj  (5  U) 

This  leads  to  an  interesting  result;  the  variance  of  the 
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phase  is  the  inverse  of  the  SNR: 


(3.13) 


Translating  this  result  to  an  uncertainty  relationship  the 
phase  uncertainty  may  be  expressed: 

AW  fct  *  iflEF  (3.W) 

'  A 

The  uncertainty  in  phase  is  the  product  of  the  time  and 
frequency  uncertainties.  This  leads  to  an  expression  for 
the  range  resolution  as  a  function  of  the  SNR. 

(3.16) 

Wide  band  ranging  systems  measure  range  by  calculating 
the  propagation  delay  of  the  signal.  In  one  such  system  a 
high  speed  code  is  transmitted.  The  received  signed,  is 
correlated  with  the  original  coded  signal  ih  a  delay  locked 
loop.  The  value  of  the  control  signal  in  the  loop  is 
proportional  to  the  target  range.  Obviously  the  resolution 
is  only  as  good  as  the  the  period  of  one  of  the  code 
bits (chips).  Other  wide  band  systems  use  other  signals  for 
ranging (chirps ,walsh  functions)  to  obtain  high  range 
resolution,  [/o]  ,  [u) 

The  system  implemented  here  at  the  Graduate  Research 
center  of  the  Moore  school  is  a  coherent  amplitude-phase 
measurement  facility.  Ranging  techniques  that  may  be 
integrated  into  this  system  are  therefore  of  special 
interest.  The  phase  of  a  point  scatterer  is  a  linear 
function  the  frequency.  This  directly  corresponds  to  the 


phase  factor  preceeding  the  scattering  integral  in  eq  (  ) . 
If  the  target  consists  of  multiple  scattering  centers;  then 
each  of  these  will  be  represented  as  a  delta  function  in  the 
in  the  reconstructed  image.  This  suggests  a  technique  using 
Fourier  analysis  to  determine  range.  First  place  the 
reference  target  in  the  microwave  field  and  measure  the 
phase/amplitude  response  over  as  wide  a  frequency  range  as 
possible.  Then  inverse  transform  this  one  dimensional 
collection  of  data.  This  will  transform  to  essentially  a 
delta  function  occurring  at  the  time  corresponding  to  the 
propagation  delay.  This  will  occur  when  the  target  has  a 
single  scattering  center  such  as  a  sphere.  Even  when  the 
object  is  more  complex  the  inverse  transform  will  be 
centered  around  the  transit  time.  In  this  Fourier  technique 
for  range  determination  the  resolution  (Ax)  is  inversely 
related  to  the  frequency  sweep  width. 

Ax  =  ~Tf  (31b) 

The  factor  of  1/2  results  from  the  fact  that  the  range  is 
half  the  signal  propagation  path  length. 

A  factor  to  be  considered  is  that  different  scattering 

centers  are  visible  from  different  receiver  positions.  It 

is  imperitive  that  all  when  several  receivers  are  used 

simultaneously  that  all  choose  the  same  scattering  center  a* 

the  phase  center  for  range  removal.  If  the  various  p» 

centers  do  not  coincide  ,  the  fringes  of  the  holoqr  «r>  . 

be  skewed.  The  data  sets  from  each  of  the  receive: 

brought  into  alignment  using  an  adaptive 
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Information  is  lost  in  the  hologram  if  the  phase  centers  for 
the  scan  lines  are  separated. 


When  this  process  is  automated  the  sweep  is  done  by 
measuring  the  response  at  discrete  frequency  points.  The 
amplitude  and  phase  are  stored  at  N  frequency  points  in  the 
sweep  range.  The  range  at  which  aliasing  will  occur  is 
given  by: 


< S./7 ) 


This  system  exhibits  processing  gain  ;  a  quality  of 
all  systems  which  spread  a  baseband  signal  into  a  wide 
spectrum.  For  this  system  the  processing  gain  is  a  function 
of  the  number  of  measurements  and  the  sweep  width. 


<a«> 


As  an  example,  for  256  measurement  points,  this  would  give 
24  db  of  processing  gain. 

A  set  of  subroutines  was  written  to  test  this  range 
removal  technique.  In  one  of  them  ,  RANGE,  the  range  of  the 
strongest  scattering  center  is  calculated  and  in  the  other 
, RANCOR,  the  phase  factor  is  calculated  and  removed  from  the 
target  data.  These  subroutines  are  used  in  programs 
SPHERE , SPHER2  ,and,  SPHER3 . 

The  time  domain  equivalent  of  this  technique  is  fitting 

a  linear  trend  to  the  phase  signal  from  the  network 

analyzer.  Since  the  phase  signal  is  modulo  (2v);this  means 
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estimating  the  frequency  of  a  ramp  waveform;  either  using  a 
least  squares  approximation  or  implementing  the  equivalent 
of  a  phase  locked  loop.  The  slope  of  the  ramp  waveform  is 
proportional  to  the  range  of  the  target.  The  accuracy  to 
which  the  range  may  be  determined  depends  on  the  sweep 
width,  the  target  structure  and  the  noise  in  the  system.  If 
the  sweep  is  wide,  then  there  will  be  more  data  with  which 
to  estimate  the  slope.  Noise  in  the  data  will  obviously 
interfere  with  the  estimation  process  as  will  as  any  phase 
shifts  due  to  the  target  structure. 

Another  type  of  system  for  generating  a  reference 
signal  utilizes  a  Target  Derived  Reference.  This  system  has 
been  extensively  studied  at  the  Electro-Optics  and  Microwave 
Optics  laboratory. [rt]  A  brief  review  of  the  ideas  developed 
to  date  in  this  regard  are  given  below.  The  complex 
exponential  in  the  integral  term  of  scattered  signal  remains 
constant  when  the  target  is  small  relative  to  the 
illuminating  signal  wavelength. 

Kp)*j  (3.1 9) 


This  is  true  if  (  )  is  sufficiently  small  such 

that(  £  •  ir  ««•!).  In  this  case  the  integral  value  approaches 

a  constant  multiplied  by  a  linear  phase;  i.e.  a  point 

scatterer.  It  will  only  occur  when  the  target  dimensions 

are  less  than  a  tenth  wavelength  of  the  illuminating  signal, 

placing  it  in  the  Rayleigh  scattering  region.  The  TDR 

signal  is  mixed  harmonically  with  a  phase  locked  scattered 
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estimating  the  frequency  of  a  ramp  waveform;  either  using  a 
least  squares  approximation  or  implementing  the  equivalent 
of  a  phase  locked  loop.  The  slope  of  the  ramp  waveform  is 
proportional  to  the  range  of  the  target.  The  accuracy  to 
which  the  range  may  be  determined  depends  on  the  sweep 
width ,  the  target  structure  and  the  noise  in  the  system.  If 
the  sweep  is  wide,  then  there  will  be  more  data  with  which 
to  estimate  the  slope.  Noise  in  the  data  will  obviously 
interfere  with  the  estimation  process  as  will  as  any  phase 
shifts  due  to  the  target  structure. 

Another  type  of  system  for  generating  a  reference 
signal  utilizes  a  Target  Derived  Reference.  This  system  has 
been  extensively  studied  at  the  Electro-Optics  and  Microwave 
Optics  laboratory,  [»»]  A  brief  review  of  the  ideas  developed 
to.  date  in  this  regard  are  given  below.  The  complex 
exponential  in  the  integral  term  of  scattered  signal  remains 
constant  when  the  target  is  small  relative  to  the 
illuminating  signal  wavelength. 

W) eirT4r  *  (S.it) 


This  is  true  if  (  )  is  sufficiently  small  such 
that(  f  •*«*!).  In  this  case  the  integral  value  approaches 
a  constant  multiplied  by  a  linear  phase;  i.e.  a  point 
scatterer.  It  will  only  occur  when  the  target  dimensions 
are  less  than  a  tenth  wavelength  of  the  illuminating  signal, 
placing  it  in  the  Rayleigh  scattering  region.  The  TDR 
signal  is  mixed  harmonically  with  a  phase  locked  scattered 


signal  at  the  imaging  frequency.  here  the  two  signal 
sources  are  phase  locked  by  a  phase  synchronizer.  The  high 
frequency  sweeper  acts  as  the  slave  signal  source.  These 
two  signals  simultaneously  illuminate  the  target.  Harmonic 
mixing  of  the  suitably  limited  TDR  and  imaging  signals  will 
yield  the  desired  phase  corrected  data. 

Another  technique  which  is  useful  for  ranging  is  the 
frequency  displaced  reference  method.  Here  the  carrier  is 
displaced  (At).  It  is  important  that  the  object  not 
contribute  to  the  phase  shift;  hence  the  displacement  must 
satisfy  the  following  condition: 


hi  ** 


2JL 
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In  this  case  the  target  structure  will  not  contribute  to  the 
net  phase  shift.  The  range  to  the  target  phase  center  is 
given  by: 

R  --  (3.ZI) 

C  *  li<jVt 

Where  (*0)  is  the  change  in  phase  for  the  carrier  and 

displace  carrier  signals  respectively.  If  the  frequency 

shift  is  small  then  the  phase  shift  will  not  be  large.  The 

resolution  of  the  system  then  is  directly  related  to  the  SNR 

in  the  receiver  channel  since  the  accuracy  to  which  the 

phase  may  be  measured  is  a  function  of  the  channel  noise. 

If  the  displacement  is  large  then  the  phase  shift  is 

increased  and  the  SNR  requirements  on  the  signal  for  a 

specific  resolution  is  decreased.  This  relationship  may  be 

expressed:  ^  r+t  m  V&N.8  < 

**i«t»w*  A  '  a? 
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The  factor  of  2  comes  about  due  to  the  phase  uncertainty 

existing  in  both  the  carrier  and  displaced  carrier  signals. 

An  alternate  method  for  implementation  of  the  displaced 

frequency  ranging  system  is  a  swept  frequency  chirp  system. 

The  ranging  signal  frequency  is  given  by  s 

,  v  *13 

(to  »/.(»♦  *0  **  ~T~ 

T-  Sweep  period 

Initial  frequency  for  sweep 
final  frequency 

The  scattered  signal  from  the  target  is  given  by: 

f,M  -■  /.'It-  (S.ZM) 

where  (dt/d^ )  is  the  change  in  frequency  due  to  the  target 
structure.  The  range  of  the  target  may  then  be  simply 
calculated  using  a  frequency  counter  and  sweep  time  T  and 
sweep  width  (f-  fj  .  This  method  could  be  used  in  an  analog 
imaging  system. 

Other  Target  Derived  Reference  systems  simulate  the  low 

frequency  reference  carrier  by  measuring  the  change  in  phase 

of  the  imaging  frequency  over  a  narrow  band.  Over  this 

small  band  the  phase  shift  is  assumed  to  be  linear.  In  one 

system  a  series  measurements  is  made  for  each  frequency 

point.  The  first  displaced  down  by  a  small  amount,  (af)» 

the  second  at  fj  and  the  final  measurement  at  (f+*f) .  Phase 

and  amplitude  are  measured  at  each  frequency  and  processed 

to  obtain  the  target  range.  In  a  similar  system  the  three 
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signals  are  transmitted  simultaneously  by  amplitude 
modulation  of  the  carrier.  These  systems  are  presently 
under  intense  investigation  by  other  workers  at  the 
Electro-Optics  and  Microwave  Optics  Laboratory  of  the  Moore 
School . 

3.2  Practical  considerations _f or  range  phase  removal 

Several  factors  influence  which  of  these  systems  would 
be  of  value  in  a  long  range  imaging  radar  system  versus  a 
controlled  laboratory  environment.  The  distribution  of  the 
reference  signal  for  complex  field  amplitude  measurement 
makes  implementation  of  system  requiring  a  central  reference 
difficult  to  implement.  Techniques  are  being  considered  in 
the  E.O.  laboaraty  for  reference  distribution  using  fiber 
optic  that  might  remove  this  limitation.  Reference 
distribution  is  accomplished  in  the  lab  readily  since  the 
distances  are  small. 

Typical  transmitter-receiver  pattern  arrangements 
might  be  the  Wells  array  [Ml],  an  orthagonal  pattern  of 
receivers  and  transmitters,  a  circular  array  of  receivers 
with  central  transmitter  or  a  random  array.  Each 
combination  of  receiver  and  transmitter  contributes  another 
line  in  the  frequency  domain  3D  data  volume.  Por  this 
reason  it  is  advantageous  that  all  combinations  of  the 
receivers  and  transmitters  are  utilised  for  data  collection. 
Reference  signal  distribution  difficulty  therefore  leaves 
the  TDR  systems  as  the  only  practical  alternatives  for  long 

range  imaging  systems.  The  AM  TDR  system  eliminates 
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reference  distribution  by  transmitting  the  reference  signal 
along  with  the  imaging  signal  and  automatically  corrects  for 
the  target  range.  TDR  systems  have  the  additional  advantage 
that  they  have  immunity  to  turbulence  and  inhomogeneities  in 
the  propagation  medium  since  both  the  reference  and  imaging 
signals  follow  the  same  path. 

There  are  several  considerations  for  determining  the  best 
TDR  system.  Narrow  band  systems  yield  only  a  weighted 
average  of  the  range  to  the  phase  center  of  the  object  while 
wide  band  system  can  resolve  individual  scatters  on  the 
target  body.  A  wide  band  system  could  adaptively  choose  one 
of  the  scattering  centers  for  the  phase  reference  of  the 
system.  A  narrow  band  system  could  not  do  this,  and  any 
error  would  introduce  image  distortion.  The  narrow  band 
systems  have  the  advantage  of  automatically  correcting  the 
target  data  for  the  range  phase  factor. 

For  the  laboratory  imaging  experiments  a  TDR  system 
would  not  yield  the  correct  phase  factor  since  in  this 
arrangement  the  object  rotates  about  an  axis  and  the  correct 
phase  factor  would  be  a  constant,  representing  the  phase 
shift  to  the  axis  of  rotation  ,not  the  target.  The  TDR 
system  looks  only  at  the  range  to  the  strongest  specular 
reflector  on  the  target  surface. 

If  movement  of  the  target  is  utilized  for  aperture 
synthesis  then  only  one  receiver-transmitter  pair  is 
required  for  3-D  imaging  and  the  adaptive  system  is  not 
required.  ,, 


Knowledge  of  the  placement  of  the  data  in  the  3-D 
frequency  domain  volume  is  necessary  for  the  reconstruction 
of  the  hologram.  The  azimuth  and  elevation  angle  of  the 
target  can  be  obtained  from  a  conventional  radar  located  at 
a  central  location  where  the  data  processing  and 
reconstruction  is  taking  place. 


IV  SYSTEM  IMPLEMETATION  OF  SWEPT  FREQUENCY  IMAGING 

This  section  o£  the  thesis  will  describe  the  research 
that  was  done  in  order  to  obtain  a  clear  understanding  of 
system  performance.  Following  this  will  be  a  section  of  the 
thesis  devoted  to  the  experimental  verification  of  the  swept 
frequency  imaging  theory.  I  ]  The  theory  is  applied  in  the 
simulation  of  the  experiments  performed. 

4.1  System  repeatability 

An  important  parameter  of  system  performance  is  the 
repeatability  of  an  experimental  measurement.  The  frequency 
range  over  which  this  possible  for  the  equipment  used 
indicates  the  bandwidth  for  which  imaging  is  possible. 
There  are  several  feedback  loops  in  the  system  which  allow 
it  to  track  variations  in  the  transmitted  signal  level.  The 
traveling  wave  tube  amplifier  characteristic  is  shown  in 
fig. 4.1.  This  plot  is  on  a  logarithmic  scale,  indicating 
that  the  TWT  amplifier  gain  drops  off  exponentially  below 
7.0  GHz  and  above  15.0  GHZ.  This  measurement  was  done  by 
first  measuring  the  system  response (cables, 
connectors, attenuator)  less  the  TWT  amplifier  and  then 
subtracting  this  response  from  the  TWT  amplifier  plus  system 
data.  A  sample  of  the  amplifier  output  is  sampled  using  a 
20  db  directional  coupler  and  this  is  fed  to  the  RF  input  of 
the  HP  8743A  reflection-transmission  (R-T)  unit.  A  crystal 
detector  at  the  'unknown'*  port  of  the  R-T  unit  rectifies  a 
portion  of  this  signal.  The  detector  output  is  brought  to 
the  external  signal  leveling  input  of  the  HP  8620C  sweeper. 


<b) 


Fig.  4.1  a)  Traveling  Wave  Tube  amplifier  gain 
character istic  in  db;  4.0-16.0  GHz.  b)  Varian  VA  618G  TWT 
amplifier  (bottom)  and  HP  8743A  ref lection-transmission  unit 
(top)  . 
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The  leveling  circuit  of  the  sweeper  can  level  the  output 
over  a  20  db  range.  In  addition  the  AGC  in  the  HP  8410B 
network  analyzer  can  track  the  reference  signal  amplitude 
over  a  40  db  dynamic  range.  System  performance  may  be  seen 
in  fig. 4. 2  for  a  cylindrical  target  7  meters  distant  from 
the  receiving  and  transmitting  antennas.  Two  consecutive 
measurements  of  the  target  were  made  and  the  results 
divided.  The  ideal  response  would  be  0  db  flat  amplitude 
and  0  degree  phase  difference  over  the  entire  frequency 
sweep.  With  few  exceptions  due  to  phase  noise  at  the  (+/-  ) 
transition  point,  the  system  has  the  desired  iepeatabili ty 
in  the  5.5  to  16.0  GHz  range.  Below  4.5  GHZ  there  are  phase 
errors  due  to  insufficient  reference  power  whereas  above  15 
GHz  errors  come  about  due  to  the  low  amplitude  of  the 
received  signal.  Noise  may  be  cancelled  by  taking  multiple 
measurements  and  finding  the  mean.  These  results  indicate 
the  useful  data  can  be  recorded  in  the  5.5-16.0  GHz  range 

4.2  Computer  control  of  target  rotation _ 

When  implementing  a  frequency  diversity  system  with 
just  one  pair  of  receiving  and  transmitting  antennas,  the 
target  must  then  be  rotated  in  the  electromagnetic  field  and 
the  scattering  measured  for  different  rotation  angles.  The 
target  used  in  the  experimental  system  rotates  on  a  stepper 
motor  driven  pedestal.  The  column  of  the  pedestal  is  1  1/2 
meters  in  length  and  is  made  of  styrofoam  material  with 
minimal  cross  section.  A  stepper  motor  controls  table 
rotation  precisely.  In  order  for  the  pedestal  to  rotate  one 


Fig.  4.2  System  repeatability  for  two  consecutive  scattering 
measurements  of  a  80  cm  long  cylinder  7  cm  in  diameter; 
4.0-16.0  GHz. 


revolution  the  motor  must  be  stepped  10  000  times.  The 
motor  is  under  direct  computer  control  using  the  digital 
output  port  the  MINC.  A  FORTRAN  subroutine  STEP 2  was 
written  for  control  of  the  stepper  motor.  It  calculates  the 
number  of  steps  required  for  a  specified  angular  rotation 
and  moves  the  table  clockwise  or  counter  clockwise  based  on 
the  direction  parameters  passed  in  the  subroutine  call.  The 
stepper  and  pedestal  are  shown  in  fig. 4. 3. 


4.3  Sphere  Simulation 

For  the  calibration  of  a  radar  system  a  reference 
target  is  required.  The  most  commonly  used  reference  target 
is  the  conducting  sphere  since  its  high  degree  of  symmetry 
does  not  favor  any  particular  polarization  for  the  incident 
illumination.  Both  the  bistatic  and  monostatic  scattering 
of  a  metallic  sphere  was  simulated. 

The  general  solution  for  the  plane  wave 
electro-magnetic  scattering  of  the  sphere  was  first  done  by 
Mie  in  1908.  [  is ]  ,  [  it]  In  the  far  field  approximation,  the 
scattered  field  is  given  by: 

[cos*  V*)©  J  M 
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S,  (e)  and  (©)  are  called  the  complex  far  field 
amplitudes  for  the  ( £ )  and  ( £ )  polarizations  respectively. 
The  quantity  in  the  square  brackets  of  eq.  4.1  is  called 


the  scattering  function. 
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The  scattering  cross  section  in  any  arbitrary  polariztion 
(  ^  )  for  an  incident  wave  polarized  in  the  (  ^  )  direction 
may  be  written: 

-  -3?- 1  (is) 


Where  (  q  )  is  the  polariztion  of  the  incident  wave  and 
( r  )  is  the  polariztion  vector  of  the  receiving  system. 
For  the  perfectly  conducting  sphere  the  coefficients  A^and  B, 
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j,(k#a)  -  Spherical  Bessel  function 
h*(k,a)  -  Spherical  Hankel  function 
( x )  -  Associated  Legendre  function 

k0  -  Wave  number  of  incident  wave 

a  -  Sphere  radius 

The  prime  on  the  expression  for  BN  denotes  differentiation 
with  respect  to  (k.  a)  . 

Polynomial  approximations  exist  for  the  Mie  series 
exact  solution. [  is  ]  Different  polynomials  are  used  for  the 
three  frequency  regions  for  scattering.  These  are:  low 
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frequency  or  Rayleigh  region  (k.a)<,4  ;  the  resonance  region 
.4<(k.  a)<20  ;  and  finally  the  high  frequency  or  physical 
optics  region  (k.a)>20.  Two  programs  BISCAT  and  SPSCAT 
implement  both  bi static  and  monostatic  cases  in  the  three 
frequency  regions.  Fig. 4. 4  shows  the  monostatic  scattering 
of  the  sphere  as  calculated.  This  is  exactly  the  same 
answer  for  the  scattering  of  the  sphere  as  the  exact 
solution.  The  horizontal  axis  is  in  terms  of  the 
dimensionless  quantity  (k6a).  Figure  4.5  a,brc  and  d  show 
the  bistatic  scattering  of  the  metallic  sphere  of  bistatic 
angles  of  30*,  60*,  90*,  and  12(?  degrees.  Note  that  the 
approximation  are  only  valid  in  the  range: 

$  <•  m-l 

where 

J 

which  leads  to  the  discontinuities  for  small  (k  a)  at  large 
bistatic  angles.  BASIC  programs  BIDISP  and  SDISP  generate 
the  graphs  for  the  sphere  simulations.  These  programs  are 
in  appendix  II. 

An  important  result  from  these  simulations  is  that  at 
high  frequencies  the  scattered  signal  is  of  constant 
amplitude  and  linear  phase  irrespective  of  the  bistatic 
scattering  angle.  The  only  exception  to  this  is  the  forward 
scattering  case  when  (©)  equals  ( V)  ,  where  the  cross 
section  grows  without  bound  as  k  increases.  This  indicates 
that  the  only  portion  of  the  sphere  that  is  scattering  for 
large  (k#a)  is  the  front  face  closest  to  both  the  receiver 
and  transmitter;  and  therefore  a  ray  optics  approximation 


Fig.  4.4  Monostatic  scattering  for  the  perfectly  conducting 
sphere.  (k  a)  varies  from  .2  to  10.5  which  corresponds  to 
the  scattering  of  a  20  cm.  diameter  sphere  in  the  frequency 
range  of  .1  to  5.0  GHz.  Log  normalized  cross 
section:  log(  ). 
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Fig.  4.5  Bistatic  scattering  of  a  20  cm.  diameter 
conducting  sphere  in  the  frequency  range  .1  to  6.0  GHz; 
.2<(k#  a)<12.6;  polarization  of  the  receiver  equal  to 
scattered  wave  polarization.  a)  Geometry  for  scattering 
expression.  b)  Bistatic  angle  30?  c)  Bistatic  angle 

60^  d)  Bistatic  angle  89*.  e)  Bistatic  angle  120? 


Fig.  4.5  (contd.)  Bistatic  scattering  of  a  20  cm.  diameter 
conducting  sphere  in  the  freequency  range  .1  to  6.0  GHz; 
.2< (k*a)<12.6;  polarization  of  the  receiver  equal  to  the 
scattered  wave  polariztion.  a)  Geometry  for  scattering 
expression.  b)  Bistatic  angle  30*.  c)  Bistatic  angle 

60.  d)  Bistatic  angle  89*.  e)  Bistatic  angle  120! 
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Fig.  4.5  (contd.)  Bistatic  scattering  of  a  20  cm.  diameter 
conducting  sphere  in  the  freequency  range  .1  to  6.0  GHz* 
,2< (k,a) <12.6;  polarization  of  the  receiver  equal  to  the 
scattered  wave  polariztion.  a)  Geometry  for  scattering 
expression.  b)  Bistatic  angle  30*.  c)  Bistatic  angle 

60.  '  d)  Bistatic  angle  89*.  e)  Bistatic  angle  120! 


may  be  applied  to  find  the  scattered  field. 
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for  the  bistatic  case.  In  the  monostatic  case  this  reduces 
to: 


Ro) 


(4?) 


For  targets  with  features  larger  than  a  few  wavelengths  in 
size  resonance  effects  become  minimal. 

Another  possible  reference  target  is  the  long  cylinder 
(1  >>  a) .  The  scattering  of  the  cylinder  in  the  high 
frequency  region  when  it  is  oriented  vertically  yields  an 
answer  similar  to  that  of  the  sphere.  For  (k  a) >5  where  a 
is  the  cylinder  radius,  resonance  effects  disappear  and  the 
copolarized  scattered  field  is  given  by: [  I S  ] 

r  *(^4^ 

The  equation  for  the  scattering  of  the  cylinder  is  used 
for  the  computer  simulations  of  frequency  swept  holography. 

4.4  Simulation  of  frequency  swept  imaging 

A  series  of  frequency  swept  hologram  simulations  were 
done  of  targets  that  would  later  be  imaged  experimentally. 
The  basic  arrangement  consists  of  separate  receiving  and 
transmitting  antennas  which  measure  the  scattering  of  a 
target  that  rotates  about  an  axis.  The  center  of  rotation 
is  chosen  as  the  phase  center  of  the  imaging  system.  For 
this  configuration  the  frequency  domain  data  lies  in  a  plane 


perpendicular  to  the  axis  of  rotation.  Therefore  the 
transforms  of  the  holograms  will  be  slices  in  this  plane. 
The  first  object  hologram  simulated  was  comprised  of  two 
cylinders  equidistant  from  the  rotational  axis.  This  target 
is  shown  in  fig. 4.6.  Approximations  for  the  various 
distances  were  derived: 

r,  2:  r-g  (&A. -e)  (H.  II*)  r  ♦  r^(|-e)  friib) 

2  r  +  *■  ©)  '  (**.HC^  r-  +e)  (l.ftd) 


the  waves  striking  cylinders  Cl  and  C2  are  given  by: 
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The  scattered  waves  from  the  two  cylinders  including  the 
cylinder  response  then  follows: 

*St i  *  Ec,  e  i A.(r,-2ft.cos|.) 

This  is  further  simplified  by  combining  terms  : 

fyp  £c©i  ^04  t-os  \  S\*&j  (If  Jlf  y 

Finally  take  the  real  part  of  this  function  for  display: 

*  C-  &•«  °',n *)  (*l‘IHb) 


This  was  done  for  a  two  cylinder  target  with  cylinders  5  cm. 
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in  radius  , separated  by  25  cm.  using  program  CYLIN.,  The 
results  were  displayed  by  program  CDISP  on  a  Tektronix  606A 
CRT  display.  CDISP  gives  the  option  of  varying  the  gray 
scale  compression  of  the  hologram  either  logarithmically  or 
by  constant  multiplication.  These  programs  are  listed  in 
appendix  II.  The  resultant  hologram  and  the  reconstructions 
obtained  through  Fourier  transformation  on  the  optical  bench 
appear  in  figs. 4.7  a,b  The  hologram  simulated  a  sweep 
from  2.0  to  18.0  GHz  in  64  frequency  steps.  The  target  in 
the  simulation  rotated  360  degrees  in  128  steps. 

Another  target  simulated  which  did  not  have  the  symmetry 
of  the  first  target  was  comprised  of  two  cylinders  both 
mounted  to  one  side  of  the  rotational  axis,  as  shown  in 
fig. 4. 8.  Two  simulations  were  done  of  this  target  with 
varying  diameter  cylinders.  In  the  first  case  7  cm  radius 
cylinders  were  used.  The  hologram  for  this  case  and  the 
Fourier  transform  reconstructions  are  shown  in  figs. 4. 8 
b,c,d.  For  the  second  simulation  the  target  was  two 
cylinders  3.5  cm  in  radius.  In  both  cases  the  cylinders 
were  located  10  cm  from  the  center  of  rotation  and  the 


simulation  was  for  a  2.0  to  18.0  GHz  sweep.  The  hologram 
and  the  transformed  images  are  shown  in  figs. 4. 9  a,b,c.  The 
two  cylinder  off  axis  target  was  simulated  by  first 
calculating  the  copolarized  scattered  field  for  a  single 


cylinder: 
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Pig.  4.7  a)  Simulation  Hologram  of  two  cylinder  target;  10 
cm.  diameter,  25  cm.  apart.  Frequency  range:  2.0-18.0 
GHa;  128  lines,  64  points/line.  b)  Optical  Fourier 
transform  of  hologram. 
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Fig.  4.8  a)  Geometry  of  two  cylinder  off-axis  target,  b) 
Hologram  simulation  two  cylinder  off-axis  target;  2.0-18.0 
GHz;  64  points/line;  128  lines,  b)  Transform  with  zero 
order  term,  c)  Transform  with  zero  order  removed. 


Fig.  4.9  a)  Hologram  simulation  of  off-axis  two  cylinder 
target  with  cylinders  7  cm.  in  diameter;  2.0-18.0  GHz  b) 
Optical  Fourier  transform  with  zero  order.  c)  Transform 
with  zero  order  term  removed. 


taking  the  real  part: 

R{^}  =  2E.  wsV^we  ♦  f &a) 

:  IE,  cos  (-\za.  — t^J s U*. a  ) 

For  the  two  cylinder  off  axis  target,  one  cylinder  is  at 
(@)=tf  and  the  other  at  (e)=90*,  therefore  the  scattered 
field  is  given  by: 

R{E4J-  c  [  cos  ( A.  *«*>©}}  fi.lfe) 

In  general  for  an  arbitrary  set  of  circular  scatterers  of 
radius  a  and  distance  1  from  the  origin;  the  scattered  field 
may  be  written: 

RfEj  ■!  Cos  {A  2  (H.n) 

This  gives  the  capability  to  simulate  the  scattering  of  any 
target  given  that  it  can  be  decomposed  into  N  spherical 
scattering  centers. 

4 . 6  Exper imental  results 

An  experimental  system  for  the  implementation  of  swept 

frequency  imaging  was  setup  in  the  anechoic  chamber  at  the 

Graduate  Research  Center  in  the  Moore  School.  The  frequency 

range  for  these  experiments  was  from  6.3  to  16.0  GHz  in  64 

discrete  steps.  The  targets  were  rotated  360  degrees  in  128 

steps.  These  holograms  were  then  identical  in  form  to  the 

simulations  previously  done. 

The  system  for  error  correction  and  range  phase  shift 

removal  was  that  used  for  high  signal  to  noise  ratio 

signals.  The  reference  target  was  a  cylinder  positioned  so 
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that  its  front  face  was  located  on  the  axis  of  rotation  a & 
in  fig. 4. 10.  A  plot  of  system  response  is  shown  in 
fig. 4. 11.  This  data  represents  the  combined  characteristics 
of  the  antennas,  amplifier  ,  cables  and  clutter.  In 
addition  it  contains  the  linear  phase  shift  that  is  the 
range  phase  factor.  As  an  example  for  the  two  cylinder 
target  shown  in  4.12  the  raw  data,  magnitude  and  phase  is 
shown  in  fig. 4. 13.  Figure  4.14  shows  how  this  data  has  been 
corrected  for  range  phase  and  system  response.  This  data 
was  generated  using  the  Fortran  program  SPHER3 . 

The  experimental  properties  of  the  two  cylinder  target 
were  studied  extensively.  Both  the  scattering  as  a  function 
of  frequency  for  a  specific  orientation  of  the  target  and 
the  scattering  as  a  function  of  angular  rotation  at  specific 
frequencies  was  obtained.  In  figs.  4.15  a,b,c  are  shown 
the  corrected  frequency  response  of  the  targe  for 
orientations  of  45*, 90*  and  135*  degrees. 

Another  computer  program  ANTPAT  was  written  to  obtain 
the  radiation  pattern  of  an  arbitrary  target  or  antenna.  In 
the  two  cylinder  case  the  pattern  was  measured  at  5.0,10.0 
and  15.0  GHz.  Note  that  when  the  cylinders  are  collinear 
all  that  is  seen  is  the  front  surface  specular  reflection  of 
the  one  cylinder  hence  the  pattern  of  a  point  scatterer  in 
the  vicinity  of  0* degrees.  These  patterns  are  shown  in 
figs. 4. 16  a,b,c.  At  high  frequencies  the  lobe  spacing  is 
much  closer  than  at  low  frequencies,  consistent  with  the 
theoretical  result  for  the  pattern.  To  see  this  examine  the 
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Fig.  4.10  Reference  target  on  pedestal?  80  cm.  long 
cylinder,  7  cm.  in  diameter. 


Fig.  4.11  Reference  target  response  including  system 
response  and  range  phase  shift;  6.3-16.0  GHz. 


Pig.  4.14  Corrected  two  cylinder  target  data  using  system 
response  o£  Fig.  4.11. 


Pig.  4.15  a)  Corrected  two  cylinder  symmetrical  target 
response;  6.3-16.0  GHz,  (e)=*yj\  b)  (e)  =  90*.  c) 
(  «  )*  I3S°. 


4.16  Scattering  pattern  of  two  cylinder  target  of  Pig. 
4.12  at  different  frequencies  as  a  function  target  rotation 
angle.  a)  5.0  GHz.  b)  10.0  GHz.  c)  15.0  GHZ. 


expression  for  the  raonostatic  scattering  of  the  symmetrical 
two  cylinder  target: 

Ej-  C  JV**1*"  .  cw-M-  (U9) 

The  expression  for  the  scattering  of  the  target  at  a  given 
frequency  as  a  function  of  angular  rotation  may  be  written: 

As  for  the  swept  frequency  response;  1  sin(e)  remains 
constant  and  the  hence  the  swept  frequency  response  is 
sinusoidal  with  period  dependent  on 

(  i( )  and  1 . 

The  final  test  for  the  system  was  the  generation  of 

actual  holograms.  The  first  target  measured  was  the  two 

cylinder  target  shown  in  fig. 4. 17a.  The  cylinders  are  of 

aluminum,  80  cm.  in  length  and  7  cm.  in  diameter.  The 

real  part  of  the  corrected  swept  frequency  data  in  the  range 

6.3  to  16.0  GHz  was  stored  and  displayed  on  the  CRT.  The 

targets  were  rotated  360  degrees  in  128  steps  yielding  a 

total  of  8192  points  in  the  hologram  (64  points/line  *  128 

lines) .  The  center  of  the  hologram  is  at  0  HZ  with  radial 

distance  directly  proportional  to  frequency.  An  example  is 

shown  in  fig. 4. 17b  and  the  reconstructions  in  figs. 4. 17c  and 

d.  These  Fourier  transforms  where  done  optically  .  [|7] 

This  procedure  was  followed  for  other  targets  not 

having  the  symmetry  of  the  first  object  used.  The  target 

type  was  the  same  as  the  simulations  done  previously.  The 

47 


‘  :t 1 


Fig.  4.17  a)  Two  cylinder  target  in  anechoic  chamber  on 
rotating  pedestal.  b)  Hologram  of  target  measured  between 
5.3  and  16  GHz  corrected  for  range  and  system  response;  128 
lines,  64  points/line,  c)  Optical  Fourier  transform  of 
hologram,  c)  Optical  Fourier  transform  without  zero  order 
term. 


first  of  these  was  a  single  cylinder  mounted  off  axis  as 
shown  in  fig. 4. 18a.  This  cylinder  was  the  same  as  the 
others  used  and  the  frequency  range  and  angular  sweep  were 
identical  to  that  of  the  two  cylinder  target.  The  hologram 
and  reconstructions  are  shown  in  figs. 4. 18  b,c,d.  The  final 
target  was  the  two  off  axis  cylinder  target  pictured  in 
fig. 4. 19a.  The  frequency  diversity  hologram  and  transforms 
are  in  figs. 4. 19  b,c,d. 


*  *9 

1 


10  Cr*\ 


Pig.  4.18  a)  Single  cylinder  off-axis  target  and  position  in 
anechoic  chamber.  b)  Experimental  hologram  of  target; 
6*3-16  GHz;  128  lines,  64  points/line;  corrected  for  range 
and  system  response.  c)  Optical  Fourier  transform  of 
hologram,  d)  Optical  Fourier  transform  without  zero  order 
term. 


7o 


i 


Fig.  4.19  a)  Two  cylinder  off-axis  target  geometry  and 
position  in  anechoic  chamber.  corrected  for  range  and 
system  response,  c)  Optical  Fourier  transform  of  hologram, 
d)  Optical  Fourier  transform  without  zero  order  term. 
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V  Conclusion 


This  thesis  has  described  an  automated  swept  frequency 
measuring  system. This  system  may  be  used  for  radar  cross 
section  measurement,  antenna  pattern  measurement  and  swept 
frequency  holography.  The  system  has  a  useful  range  of 
5.0-17.0  GHz  in  which  amplitude  and  phase  of  the  scattered 
microwaves  from  targets  in  the  anechoic  chamber  of  the 
Graduate  Research  Center  may  be  recorded  and  stored.  A  DEC 
MINC  LSI-11/2  completely  automates  the  data  acquisition 
process.  A  complete  error  correction  algorithm  was 
implemented  using  the  data  storage  and  processing 
capabilities  of  the  minicomputer. 

The  effect  of  range  phase  shift  on  swept  frequency 
holograms  was  investigated  and  various  techniques  for  its 
removal  were  investigated.  It  is  believed  that  a  TDR  system 
for  range  phase  removal  is  required  for  implementation  of  a 
practical  radar  system.  This  system  must  have  extremely 
high  resolution  for  coherent  imaging.  The  relationship 
between  TDR  system  bandwidth,  receiver  channel  bandwidth  and 
resolution  was  derived: 

-  "is- .  = 

M  sir  *1  nr  A  af 

where  (&f)  is  the  imaging  bandwidth,  (R)  the  range 
uncertainty,  B  the  receiver  channel  bandwidth  and  N.  the 
noise  power  spectral  density. 

Simulations  were  performed  for  the  scattering  of 

various  radar  targets  which  include  the  conducting  sphere, 
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infinite  cylinder  and  combinations  of  these.  The  holograms 
recorded  from  these  analytical  results  reconstruct  the 
targets  extremely  well  and  set  a  goal  for  practical  system 
performance.  An  expression  was  derived  for  the  scattering 
of  N  spherical/cylindrical  cylinders  in  a  plane  passing 
through  their  centers: 

X  ^ {***■*« (et9n)J  £4*17) 

where  1  is  the  distance  from  the  axis  of  rotation  (  en)  the 
angle  relative  to  some  reference  for  the  target  angular 
position  and  a  the  target  radius. 

Finally  experimental  swept  frequency  holograms  were 
generated  using  a  rotating  pedestal  under  computer  control 
to  scan  the  target  in  one  dimension.  The  experiments  done 
indicate  the  feasibility  of  implementing  a  practical 
holographic  radar  system.  The  holograms  obtained  for 
various  targets  agree  well  with  theory  even  though  the  error 
correction  and  range  phase  shift  removal  techniques  used 
were  robust  in  nature.  It  is  believed  that  better  images 
are  possible  given  that  the  error  correction  techniques 
previously  outlined  are  implemented. 

Further  work  can  be  done  in  testing  the  TDR  techniques 
for  their  suitability  for  an  imaging  system.  The  system  may 
also  be  expanded  to  include  scanning  in  the  (  )  direction 
to  give  true  3-D  imaging  capability.  This  may  be 
implemented  by  adding  a  stepper  motor  controlled  azimuthal 

scanner  to  the  top  of  the  rotating  column. 
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